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In these lectures we consider some topics of Quantum Field Theory in Curved Space. In the first one particle 
creation in curved space is studied from a mathematical point of view, especially, particle production at a given 
time using the so called "instantaneous diagonalization method". As a first application we study particle pro- 
duction in a no-oscillating model where re-heating may be explained from gravitational particle creation. In 
the second one we re-calculate, with all the mathematical details, particle production in the Starobinsky model. 
Particle production by strong electromagnetic fields (Schwinger's effect) and particle production by moving 
mirrors simulating black hole collapse are also studied. In the second lecture we calculate the re-normalized 
two-point function using the adiabatic regularization. The conformally and minimally coupled cases are con- 
sidered for a scalar massive and massless field. We reproduce previous results in a rigorous mathematical form 
and clarify some empirical approximations and bounds. The re-normalized stress tensor is also calculated in 
several situations. Finally, in last lecture quantum correction due to a massless fields conformally coupled 
with gravity are considered in order to study the avoidance of singularities that appear in the flat Friedmann- 
Robertson- Walker (FRW) model. It is assumed that the universe contains a barotropic perfect fluid with state 
equation p — top (being p the energy density and p the pressure). The dynamics of the model is studied for all 
values of the parameter uj, and also for all values of the two parameters, that we will call a and /?, provided by the 
quantum corrections. We will see that only the case a > could avoid the singularities. Then when uj > —1, 
in order to obtain an expanding Friedmann universe at late times (only a one-parameter family of solutions, no a 
general solution, has this behavior at late times), the initial conditions of the no-singular solutions at early times 
must be very fine tuned. These no-singular solutions are: a general solution (a two-parameter family) leaving 
the contracting de Sitter phase, and a one-parameter family leaving the contracting Friedmann stage. On the 
other hand for uj < — 1 (phantom field), the problem of the avoidance of singularities is more involved because 
if one considers an expanding Friedmann stage at early times, then instead of fine tune the initial conditions 
one also has to fine tune the parameters a and ft to obtain a behavior without future singularities, because only 
a one-parameter family of solutions follows a contracting Friedmann phase at late times, and only a particular 
solution behaves like a contracting de Sitter universe. The rest of solutions have future singularities. 
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I. INTRODUCTION 



In these lectures we'll try to give a self-consistent presentation of the quantum field theory in curved space and also 
to show some of its applications. In the first lecture we will give a mathematical presentation of the subject and we 
will re-derive, with all the details, some of its applications to the theory of gravitational re-heating. Our presentation is 
no-standard, in the sense that we start explaining the "adiabatic vacuum prescription" introduced by Parker in his thesis 
iQJ-d], and the "instantaneous diagonalization method", introduced in Russian literature at the beginning of 70's, based 
in the idea that the number of created particles in a given mode at a given time is the energy of the mode at this time 
divided by the energy of a single particle in that mode 0-01 ■ As an instructive example we calculate particle creation 
in the flat Friedmann-Robertson- Walker (FRW) chart of the de Sitter space where one can sees the difference between 
both prescriptions. After that, we introduce the "in" and "out" states in asymptotically flat spaces where particle creation 
can be defined in the standard way As an application, we re-calculate the gravitational particle production in 

a transition from the de Sitter phase to the radiation-dominated one 0], and we also discuss the problem of a second 
inflationary stage related with the back-reaction IU5I417I1 . This example is important because it describes approximately 
the inflationary phase followed by a transition to a radiation dominated universe, and the particle production process may 
be used to explain the pre-heating in inflationary non-oscillatory models ITil [Hfl, [l9tl . Finally we study, in a great detail, 
particle production in the Starobinsky model [20-22] because we believe that there isn't a clear explanation of particle 
creation in this model. In this case, in order to obtain the well-known results, firstly one must disregards the power-law 
expansion of the universe and only retains the oscillating behavior of the scale factor at late times, secondly one has to 
assume that the energy density of the created particles is a well-defined quantity in this model and then one has to choose 
a particular form of it, and finally one has to assume that one kind of particles, named scalarons, are the responsible for 
the late time behavior of the universe, and also that is the decay of these scalarons what produce particles. 

Particle production by strong electromagnetic fields l23ll are also studied, and Schwinger's formula j24Tl that calculates 
the probability that the vacuum state remains unchanged in the presence of a constant electric field is deduced, in an 
elementary but not at all mathematically correct way, using standard methods of Quantum Field Theory. More precisely, 
it is deduced calculating the Bogoliubov coefficients for every mode via the W.K.B. method in the complex plane. At the 
end a rigorous demonstration is outlined. 

The last part of this first lecture is devoted to the study of particle production by moving mirrors. Our interest is 
concentrated in trajectories that simulate the black body collapse, and to get the same kind of results obtained by Hawking 
in l25ll . i.e., to obtain the black body spectrum. This occurs for perfect reflecting mirrors, but we'll show that for semi- 
transparent moving mirrors the radiation spectrum is a bit different. 

In the second lecture the vacuum quantum fluctuations are studied. We calculate the re-normalized two-point function 
subtracting adiabatic modes up to order two. We do the calculation for conformally coupled fields and also for minimally 
coupled ones re-obtaining, in a consistent way, all the early well-known results. Calculating the two-point function is very 
important in the context of inflation, for example, the eternal inflation phenomenon is manifested in some inflationary 
model (new inflationary universe l26l - l28ll . chaotic inflation ll29ll ). i.e., the large-scale quantum fluctuations of the inflaton 
field termed by its two-point function lead to a process of infinitely self -reproducing inflationary mini-universes [30]. 
Studying the back-reaction of particles produced in the pre-heating phase also requires the two-point function of the 
inflaton field and the two-point function of the light particles involved in such a process ll3ll - [33ll . 

The mean problem with the two-point function is that it is ultra-violet divergent and requires re-normalization. The 
simplest method for obtaining divergence-free expressions is the adiabatic regularization based on subtracting some gen- 
eralized WKB modes [341 that have the same behavior at large frequencies as the exact modes; the divergent terms then 
cancel. But the procedures for calculating the re-normalized two-point function differed somewhat in early work. The 
authors assumed some not fully justified frequency cut-off or made unjustified approximations in order to obtain finite 
quantities HSHHH- 

Our task in this lecture is to matematically clarify the features appearing in those works. Firstly, the massless case is 
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studied (conformally coupled and minimally coupled case), reproducing all the previous results in full detail. After this, 
we study a massive field in the de Sitter phase where it's assumed that its mass is smaller than the Hubble parameter (this 
is typical in the inflationary models l30lo . Here we derive the two-point function at late time in full detail and accurately 
demonstrate the mean formula obtained in [26]. We finish this lecture reviewing some important results about the stress 
tensor re-nomalization which will be used in last lecture. 

In last lecture we study the avoidance of cosmological singularities if one takes into account the vacuum corrections 
due to a massless conformally coupled field. 

It's well-known that the classical solutions of the general relativity for a Friedmann-Robertson- Walker (FRW) model 
contain, in general, singularities (Big Bang, Big Rip, future sudden singularities), this means that near these singularities 
the space-time curvature is arbitrarily large. Then, for curvatures on the order of the Planck length, quantum effects 
have to be taken into account. These quantum effects, can violate the so-called energy conditions 13811 . and consequently 
they can modify drastically the classical solution. For this reason, it is possible that quantum effects avoid the classical 
singularities [39, 40]. 

We consider the quantum effects produced by massless fields conformally coupled with gravity. This is an special 
case where, for a flat FRW universe, the quantum vacuum stress tensor, that depends on two regularization parameters, 
that we call a and /3, can be calculated explicitly. Then due to the trace anomaly and the equation of conservation, one 
easily calculates the vacuum energy density that contributes to the modified Friedmann equation. This equation cannot be 
analytically integrated, but a qualitative phase-space study can be performed. This is the main objective of this lecture. 

First, we introduce the quantum effects and write the modified Friedmann equation that depends on the parameters a 
and f3, which we'll assume that can take all possible values. After that, we study the simplest case, i.e. a = 0, in this 
case the modified Friedmann equation becomes a first order differential equation and can be integrated. Our conclusion, 
in that case, is that the singularities are not avoided. Another simple case corresponds to the case of an empty universe 
(it doesn't contain any barotropic fluid, only quantum effects are taken into account). An special case (a < 0, /3 < 0) 
is the Starobinsky model 12011 . We'll see that, in that case, all solutions contain singularities, except when j3 < 0, where 
it appears an unstable de Sitter solution, and an unstable solution that connects the de Sitter solution with the point 
H = H = (being H the Hubble parameter). Finally, we study the general case, i.e., a universe filled by a barotropic 
perfect fluid with state equation p = ojp. The only case where no-singular solutions may appear is when a > and j3 < 0. 
Then taking the same point of view as j4lll42Tl . we show that when lj > —1, the no-singular early time behaviors that 
can lead, at late times, to the Friedmann expanding stage, are a contracting de Sitter phase and a contracting Friedmann 
phase. However their initial conditions can be very fine tuned in order to match with the expanding Friedmann stage. On 
the other hand, for lu < — 1 the late time behavior of no-singular solutions that come from the expanding Friedmann stage 
at early time, are the contracting Friedmann phase and the contracting de Sitter one. In this case, instead of fine tune the 
initial conditions, one also has to fine tune the parameters a and f3 in order to obtain no-singular solutions. But in both 
cases, the no-singular solutions are unstable in the sense that an small perturbation leads them to a singular behavior. 

The units used in these lectures are c = h = 1. 



H. PARTICLE CREATION BY CLASSICAL FIELDS 

Particle creation by gravitational fields is studied in this section. The theory developed is applied to a non-oscillating 
inflationary model and to the Starobinsky one. 

A. Graviational particle production 

1. Quantum fields in curved space-time: General Theory 

It's well known that the Lagrangian density of a scalar field is Hill C — ^(d^d^rfi — m 2 4> 2 — £i?</> 2 ), and its corres- 
ponding Klein-Gordon equation is given by 

(-V M V + m 2 + = 0, (1) 

where £ is the coupling constant and R is the scalar curvature. If one considers the flat FRW metric ds 2 = 
—dt 2 + C(i)dx 2 — C(rj)(—dri 2 + dx 2 ) (being t] the conformal time), the modes of the form 0k(x, 77) = 
(2-7r) _3 / 2 C _1 / 2 (?7)e lkx Xk(f7) will satisfy the equation 

Xk(»7) + ^(»7)Xk(t»)=0 J (2) 
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where we have introduced the notation = w k( 7 /) + (£ — l/®)C{r])R(r]), with w^(ry) = m 2 C(rj) + |k| 2 and 

flfa) = 3 - §g£) = 6^ (being C = a 2 ). 

Firstly, we are interested in the conformally coupled case, i.e., in the case £ = 1/6, where equation (O becomes 

Xkfa)+Wk(»?kkfa) = 0. (3) 

This is the equation of a set of no-interacting harmonic oscillators, so to built a quantum theory, we can consider the 
following Hamiltonian operator H(rf) = | (f) 2 + cj 2 (r;)0 2 ^ — corresponding to a single harmonic oscillator. In 

Heisenberg picture the operators n and 4> satisfy the equations — n and n = —uj 2 4>, that is, $ satisfy the Klein-Gordon 



equation (f>" 

where the mode function \ is a solution of the Klein-Gordon equation and A x is, in that picture, a constant operator that 
we will call the "annihilation operator relative to the mode \ "• 

From the commutation relation [fl, </>] — —i, one can deduces that x must satisfies the relation x'x* ~ x'*X = — 
which means that the annihilation operator is given by 

4 = -*(x'*(^)-x*Wn(r?))- (5) 

Once this operator has been introduced one can defines the "vacuum state relative to the mode x"> namely |0; x)> as 
the quantum state that satisfies A x |0;x) = 0. It's clear from this definition that there isn't particle production at any 
time, because (%; 0|At.A x |0; x) = all the time (43J. However, this definition depends on the choice of the mode x- 
Effectively, if one chooses two different mode functions, namely \i an d X2, since A Xl — oti t 2A X3 + /3* w ^ m 
ai t 2 = — iyV{x2', Xi] an d P12 — ~^[X2' X*] (where W denotes the Wronskian), then an observer in the |0; X2) vacuum 
state can observes xi-particles because one has A/1,2 = {X2\ 0|A^i4 Xl |0; X2) = | /?x . 2 1 2 - 

In this way, if one considers the family of solutions to the Klein-Gordon equation, namely Xrf (v)> defined by the initial 
condition 

Xri'iv') = fiv')') X^'W) = oiv')) where / and g are some arbitrary functions, (6) 

one can calculates the number of Xr-particles detected by an observer in the |0; Xt') vacuum state, that is, the number of 
produced particles at time r from the vacuum state at time r', with the formula 

M{t;t')= ( X t>;0 h \A XtjH A Xt , h \0 h ; X t>) = |/?(r;r')| 2 , with 0(t; t') = iW\Xr'\ Xt}- (7) 

Note that, different families of solutions give rise to different definitions of the vacuum state. For example, to define 
the adiabatic vacuum modes first we consider the e-Klein-Gordon equation ev" + uj 2 (j])v = 0, where here, e is a 
dimensionless parameter that one shall set e = 1 at the end of the calculations. At order "n", a WKB solution of the 
Klein-Gordon equation is (see for details |44T0 : 



y 2VV„(r;eJ 



where Wq — uj and 



W n — terms until order e 2 ™ of t u> 2 — e 



2W n -! 4 W 2 _! 



(9) 



Once one has introduced the WKB solutions, the adiabatic vacuum at order n is defined through the family Xnrf (v) 
that satisfy the initial condition 

Xn; V , {v')=Xn;WKB{r]'\^ = i y, Xn; V ' iv' ) = Xn;W K B W ^ = ^) ■ (10) 
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From this definition, at the n order, the /3-Bogoliubov coefficient is given by /3„(t;t') = iW[Xn-,T'\ Xn;r]> an d the 
number of produced particles at time r can be calculated from the formula 

Af n (r;r')^\P n (r;r')\ 2 = Q[|x^(r)| 2 + ^ 2 (t)| X „ ; .<(t)| 2 ] \w n {r) 

+ S^ IX " ;r ' (T)|2 + m^r) {x '^ AT)Xn -' T ' {T) + ^A-r)xl;Ar)), dD 

where W„(t) = W n (r; e = 1). 

Another important family of solutions to the Klein-Gordon equation, namely Xdiag;r]'(v)> is given by the initial condi- 
tion 

Xdiag- V '(v') = XQ;WKB(ji\t = 1); Xdiag: V ' W) = -i^(v')X0;WKB (v'', £ = l )- ( 12 ) 

This family defines the so-called instantaneous Hamiltonian diagonalization method, and their Bogoliubov coeffi- 
cients 

^dtagir; T') = -iyV[Xdiag-T'\Xdiag-X Pdtag{T] t') = i\V[Xdiag;T> 5 Xdiag;r], (13) 

can be calculated as follows: Writing Xdiag-y (v) = ^diagi^r^Xdiag^iv) + 0dia g {r;T')x diag . T (r])' at V = t one gets 
the system 



Xdiag-T'{T) = a d iag{T\T')xO;WKB{T\e = 1) + Attag (t; t')x5 ;VVKB (V; £ = 1) 

x'diag;A T ) = "M 1 ") (adiag(r; t')xo-,wkb(t; e = 1) - /3 dia9 (r; t')xo.wkb( t > 6 = *)) - 
which can be used to obtain the interesting formula 

Af diag (r;r') = \l3 dlag (r;r')\ 2 = Q[|x^(r)| 2 + W 2 (r)| Xdws;T ,(r)| 2 ] - i W (r)) , 



(14) 



(15) 



which shows that this quantity is the energy at time r of the mode Xdiag-.r' divided by the energy, at time r, of a single 
particle. 

Another way to obtain the Bogoliubov coefficients can be done from the system ( TBl i if one takes into account that the 
mode function Xdiag-^'ij) satisfy the equation Xdiag-T'( T ) + ^ ''( T )Xdiag-T' (t) = 0. As a function of the variable r, one 
has 

f a diag( T > T ') = u'ir) (v*. wkb (t;£ = 1)) 2 ^agi^T') ^ 

\ P'diag( T ^ T ') = w'(t) K ra (^;£ = l)) 2 a<i 1 « ! ,(T;r'). 

This system can be solved by iteration, for example, if in the first iteration one chooses adiag ( T > T ' ) — 1> one will arrive 
at formula 

Pdiag(r;r') ~ £ |M e -«/^(V)^ d?? . (17) 



Example II.l. (Particle creation in the flat FRW chart of the de Sitter space-time [45, 46]) 

In the de Sitter phase the scalar factor is given by a(r/) = —l/(Hrj), with — oo < r/ < (being H the Hubble 



parameter), and the frequency has the form uj(rj) — w + jji^i , where ujq is a constant. We are interested in the case 
m 3> H which correspond to the adiabatic approximation (see below). An easy calculation yields 

\ im/ H 



Xo.,WKB(me = l)= * WW(^) H J (18) 

which shows, when r\ — > — oo, that 



X0;WKB(m e = !) -> \l %UJ ° V ' Xo-WKBi 7 !'^ = !) -iuJ vo-,WKB(m e = !)• (19) 
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The mode solution that satisfy the initial condition ( 1791 ) is given in terms of the Hankel functions [11] 



X(V) = CJ^HW fat,), (20) 



with fi = J\- % = im/H, C = e-^+5) ^ e W(2*0 e -if and v _ e -i^^ 
Using the asymptotic form of the Hankel functions at late times, i.e., when \tj\wo < 

x(v) * -Cjf* 
an easy calculation provides that 



e-~/*r ( i - {^y m/H - r ( i + W#) ^fY" " " 



(21) 



l/WO; -oo)| 2 = \T (1 - im/H)\ 2 e^ H . (22) 

Using at this point fcfy] \T(1 + iy)\ 2 = try/ s'mh(iry\we conclude that, when \r)\uo <C 1, the number of produced 
particles, using the diagonalization method, is given by l \4al 

H 2 

Afdiag(0;-^) = T7T -j. (23) 
16m z 

However if we use the zero order ( the other orders give the same result) adiabatic vacuum modes, the square of the 
(3-Bogoliubov coefficient will be given by 

|A,(0; -oo)| 2 = \x(0)xo;wkb(0; e=l)- x' (0)Xq-,wkb(0; e = 1)| 2 . (24) 
Inserting ( 1271 ), in this last formula, we obtain 

A/" o (0; -oo) = — |r (1 + im/H)\ 2 e -™' H = ( e 2mn ' H - l) ~* . (25) 



This is the thermal spectrum obtained in the flat FRW chart of the de Sitter space-time 4451 14 
Remark II.l. The two methods give a different result because when \ti\ujo <C 1, one has 

Xo,wkb(V: e = 1) = -i^{v)xo : WKB{V: e = 1) + ^Xo,wkb(v, e = 1) 7^ -iu(r])xo,WKB(v, e = l ), (26) 

?/iis is afwe fo the fact that lim„_>.Q- w(?y) = oo, f/W is, a? iafe fime f/iere is no? a well-defined "out" region (see below for 
a precise definition of the "out" region). 



In general (for arbitrary values of £), for a given set of modes satisfying W[XkjXk] = — * one can expands the 
quantum field, in the Heisenberg picture, as follows: 0(x, 77) = ^ k A Xk </>k(x, 77) + A Xk k (x, 77), then one can defines 
the quantum vacuum state relative to the modes 0k(x, v) = (2> 7r )~ 3 / 2 C~ 1 / 2 (r])e lkx Xk(v)> namely |0; x), which must 
satisfies A Xk |0; x) = for all values of k. However if one considers another set of modes, namely 0k (x, 77), one also 
may develops the quantum field as </>(x, 7/) = ^ k A^ k ^k( x ; v) + ^^ k ^k( x ' V)' wner e 0k = «k0k + m ^ thus, 

A Xk = a'^A Xk — /?k^x k w i m l a k| 2 — |/?k| 2 = 1. The vacuum |0; x) relative to the modes </>k(x, rj), is related with the 
other vacuum trough the relation 

|0;x)=ncxp|i^(it k J 2 || ; X ), (27) 

and the operator "number of particles in the mode k" that depends on the choice of the set of modes, for example, for the 
set k (x,77) isA/" Xk = i Xk i Xk satisfies (x; 0|A/" Xk |0; x) = 0, however one has (x; 0|A/" X1C |0; x) = IA<| 2 - 
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Once we have introduced these definitions, one can considers the family of solutions to the Klein-Gordon equation, 
namely Xk;ri'{il)< defined by the initial condition 

Xk;r,' {if) = fk{rf); xL-tf (?/) = Sk(»/)> where h and g k are some arbitrary functions, (28) 

to calculate the number density of Xr-particles per unit volume detected by an observer in the |0; Xt') vacuum state, that 
is, the number density of produced particles at time r from the vacuum state at time r'. The general formula is 



N{T;T>) = (2^F / d3k ^;0|it k;T i Xk; jO; X ^) - ^3 / d 3 k|/3 k (r;r')| 2 , 



(29) 



wifhj8k(T;T , ) = iW[xk i r';Xk i r]. 



Remark II.2. In the conformally coupled case one also can defines the energy density of produced particles at time r 
from the vacuum state at time t' as follows p{r; t') = -p^jr; 

Example II.2. (Particle production in the adiabatic approximation) 

This approximation is based in the assumption fl^ <C Q 2 ,. In the conformally coupled case this assumption becomes 
ui' k -C w k and it is always satisfied when H <C m. In that case one has 



XKdiag-Ar) = Xm-.wkb{t;€ = 1) = , /—L-e^/r' (30) 

y m(t) 

and when 0J k (r') = one can inserts this expression in formula ( 1751 ) to obtain 

\Pk,diag{T,T )\ - - Wd^r) ' K ' 

which helps us to conclude that, using the instantaneous diagonalization method, the number density of created particles 
per unit volume is given by /0/ 

, mH 2 f°° x 2 mH 2 



ant/ f/ze/r energy density is 



„ m 2 H 2 f°° x 2 m 2 H 2 



Remark II.3. Note that formula ( 1231 ) can be easily obtained from formula ( 1571 ) applied to the de Sitter phase. 

To finish we consider and asymptotically flat FRW space-time, that is, we assume that lim,,_>. = b 0o C{rf) = C±, and 
we take the following set of modes <pin k(x, 77) = (27r)~ 3 / 2 C~ 1 / 2 (r7)e zkx e , k ' when r\ — > — oo, and the the set 

<Ao«t,k(x, r\) = (27r)~ 3 / 2 C~ 1 / 2 (?7)e lkx e , +,k " when 77 -» 00 (being w± k = \/ m2 C± + M 2 )- ^ n tn i s context we can 

define the "in" and "out" vacuum states, namely |0» n ) and \0 O ut), and the "in" and "out" annihilation operators Aj„ k and 
A ou t,k. Then, the average number of produced pairs, at late times, in the k mode, is given by (0in\Jvout,k\0iTi) — |/3k| 2 , 
where N out ,k = A*. V A 

out,k is the operator number of "out" particles, and the beta Bogoliubov coefficient is obtained 
trough the relation <p lnM = a k </> out >k + M* out k . 

In that context, the number density of created particle per unit volume, at late times, is given by 

N=-^ f d 3 k\^\ 2 , (34) 
(27ra) J J 

and their energy density by l49ll 

1 



(2ira) 3 a 



d 3 k w+ , k |/i k | 2 . (35) 
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In general, it is impossible to solve the mode equation (O but one can rewrite this differential equation in an integral 
one as follows JH Ha : 

Xk (v) = + / Vk(ry') sin( W _, k (ry - T]'))Xk{v'W , (36) 

where Vk(r]) = w 2 k - f2 k (r?). 

Applying Picard's method to lowest order, i.e, replacing Xkirf) by —j — =- one obtains the following approximation 
that works very well for massless nearly conformally coupled fields 

a k £n + -^— f V k ( V )dr) [ e-^- ^v^dr, (37) 

^-,k JR A^-.k JR 

Note that, one will have to assume lim^-too Vu(r]) — if one wants well defined Bogoliubov coefficients. This always 
happens in the massless case, and using Plancherel's theorem is not difficult, in the massless case, to prove that 

N= ^~ lf f I a\ V )R 2 (v)dv. (38) 
167ra J ./d 



2. Particle production in the transition from de Sitter phase to a radiation dominated universe 
Consider the following scale factor ffl 



a(r/) = { „, \ n i f ~ (39) 
\ H(t] - rjo) - -jj^ for ry < V, 

where rjo < is the time when the sudden transition occurs. 

This example is interesting because in the radiation phase massless particles cannot be produced, so in that phase, 
the number density of created particles and their energy density are well defined quantities. Moreover, it describes 
approximately the inflationary phase followed by a transition to a radiation dominated universe, and the obtained result 
can be used to explain the reheating process of the universe after inflation in some no-oscillatory models lfl9ll . 

First, we wil consider massless nearly conformally coupled particles. In the de Sitter phase, one has il k (?/) = |k| 2 + 
(12£ — 2)/r/ 2 , and in the radiation one fi k (r/) = |k| 2 . Using formula fl37T i one gets 



|k| loo V 2 



/?k — 7TT / e- 2 ^^^d V , (40) 



and form formula (l38l one easily obtains 



(6£-l) 2 = (6tj-l) 2 H 3 /any 
127ra 3 ^ 12tt \ a J 



^ = -^n^- v s ,:L (-) - (41) 



where ao = a(r/o). 

It isn't difficult to show, from formula d37| i, that the energy density diverges. However if one assumes that the transition 
is not abrupt one obtains (see 11411 ) 

P~^^^-1) 2 W^) 4 . (42) 

Now, we consider massless minimally coupled particles where the mode functions that describe the vacuum state in the 
de Sitter phase are given by 

Xk(ry) = - JZ^ff$(|k|rj) = K \— < ' k " [ I + — ) - (43) 
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and in the radiation one, these modes have the form 



(44) 



Matching at time n = tiq one obtains /3k = ol .| 3 s e 2i l k l 1 )° what implies that the number density of produced particles 

Z|K| ?? 

is infrared divergent, and their energy density is both, infrared and ultraviolet divergent. 

To eliminate the infrared divergency one can imagines that at very early times the universe is in a radiation phase, 
then at a given time, for example rj = —H^ 1 (t = 0), there is an abrupt transition to the inflationary phase described 
approximately by a de Sitter one (see for details i36H5lll ). On the other hand, to avoid the ultraviolet divergencies and 
can assumes that inflation finishes with an smooth transition to the radiation phase lfl4ll because in that case modes with 
|fc| \r]\o have a very small contribution. Anyway if one is only interested in production of particles whose modes 
leave the Hubble horizon, that is, in modes that satisfy H < Ikl < 1 770 1 1 , one can uses the formula /3k = ^, e ~ 2 'l k l'?o 
to obtain 



8tt 2 o 3 ^ J h 11 11 8tt 2 V a J \H\r) 
and 

87r 2 a 4 ^7 H 8vr 2 V a / \H\vo\J 

Now we consider the general massless case where the modes that defines the vacuum state in the de Sitter phase are 

X k (v) = cJ^H^(\k\ V ), (47) 



with v=J\ - 12£ and C = e^^ + f). 

For modes that satisfy | kr/o | <C 1 one can uses the asymptotic formula 



H^(z) = % - (z/2)-T(^) - (z/2)" — i— . (48) 



Then matching at point 77 = 770 one obtains for v > (i.e., for £ < 3/16) [5 



l/?k| 2 = Tjb (Mvo\/2\)-*'- 1 1» (1/2 - vf . (49) 



In the opposite case |kr?o| ^ 1, from the asymptotic formula |4' 



H<?\z)^C X l^(l-i^P±)e- i \ (50) 



after matching at point r\ = rjo one obtains 



7TZ V 8z 



1 \v l - 1 



16 |k| 2 r^ 



what means that the energy density is always ultraviolet divergent. 

From these results one concludes that for 5/48 < £ < 3/16 (0 < v < 1) there isn't infrared divergencies and thus the 
number density of particles converges. Moreover, if one assumes that the transition is smooth then their energy density is 
also finite. On the other hand, when £ < 5/48 (y > 1) the number density is infrared divergent (and for £ < {v > 3/2) 
their energy is also infrared divergent). The solution to the avoidance of divergencies in this last case is the same as for 
the minimally coupled case. Then if one is only interested in modes that leave the Hubble Horizon, one can uses the 
approximation |k??o| -C 1 to obtain the formulae 

4 " -2,w 1/0 ^2„ 3 /ao\V 1 



32^731) r » (1/2 -rnfj Uhol)^) - 1 '' (52) 
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and 

4" 



TJCT^T- 1 • < 53) 



' 64tt 3 (2i/-3) ww ; \aJ \(H\r]o\)^- 3 '> 

Remark II.4. Note that the case v = 1/2 « f/ze conformally coupled case and thus there isn 't particle production. 

To obtain results for massive particles we will apply the diagonalization method, more precisely, to calculate 
\Pk,diag(r]o] — o°)| 2 we have to use formula ([T5l > with, and that is very important, Wk(f?o) = ^/|k| 2 + gJj^^ ■ Then 
for light particles (particles with small mass compared with the Hubble parameter, i.e., with m <C H) no-conformally 
coupled one can uses the formulae ( l52l and (1531 1. However for conformally coupled particles we've obtained, after a 
cumbersome calculation, in the range raj H <C |k?yo| 1 

\Pk,diag(V0; "00)| 2 - g4?? 4| k |4 ' (54) 

then, integrating over the modes that leave the Hubble horizon one has 

pW%)~™ 4 (^) 4 ln(^). (55) 

In the opposite case, that is, for m 3> H, from the example 2.1 in the range |k^o | "C 1, one gets 

H 2 

IpKdiagiW, -00)| 2 ~ — (56) 

and integrating over the modes that leave the Hubble horizon one has 

Pdiagim) ~ — ( — ) ■ (57) 

Finally it's important to remark that for light particles it is well-known that its energy density during the de Sitter phase 
is of the order H 4 (see equation ( |194t in lecture II), and then if one takes a inflationary model where the energy density 
of the inflaton field, namely p v , satisfies p v <C mt (being m p the Planck mass), for example if it is at the grand unified 
theory scale 15311 . then from the Friedmann equation H 2 = one will deduce that H 4 <C p v , that is, the back-reaction 

effect will be negligible. However, when p v ~ m p , in lfl5l[l6l the authors showed, in the minimally coupled case, that 
there is a new inflationary state (which makes the results obtained above incorrect) driven by the field that produce light 
particles, namely cf>, and not by the inflaton field, namely <&, and thus instead of studying gravitational particle production, 
one should study the mechanism of production of particles of the field <f> by the oscillations of the field <f> due to the 
potential energy density ^m 2 (f> 2 . 



3. Particle production in the Starobinsky model 



In this model the universe e merg es form the de Sitter phase and at late times the scale factor is approximately 
given by (see for example ll20l l2lll ) a(t) = i 2 / 3 [l + ^^sm(Mt)] (in terms of the Hubble parameter H(t) = 



cos 2 (Mi/2) 



1 



Mt 



where M is a constant that is related with the vacuum polarization effect produced by 



massless conformally coupled particles. More precisely, this model corresponds to an empty universe (it doesn't con- 
tain any barotropic fluid), and only quantum effects due to massless conformally coupled fields are taken into account 
1I20I - I22I I4II1 . The energy density that contributes to the Friedmann equation depends on two parameters lf39L l54l l55ll . 
one of them is M and the other is called Hq in [21]. Then, when these two parameters are positive there is an unstable 
no-singular solution, that emerges from the de Sitter phase with H(t) = Hq and at late times approaches to a matter 

(A more detailed description of this model is given in 



1 



sin(Aft) 
Mt 



dominated universe with H(t) = ± cos 2 (Mi/2) 
last lecture). 

First at all we study the production of massless nearly conformally coupled particles. To do this we'll follow Vilenkin's 
viewpoint l2~lll (see also [22]). The idea is to calculate the rate of particle per unit volume and per unit time, this concept 



11 



was introduced in rf50Tl but without any explanation. Here we'll try to justify the results obtained in 1 50] . Note that to 
calculate the particle rate per unit time one needs to calculate the particle production at a given time and then to calculate 
its derivative with respect to time, that is, one has to use formula ( 1291 with t' = — oo. To perform this calculation we've 
to choose a given family of mode solutions. From our point of view the most natural choice is 

XK r(v) = ^t=S + iu r W) MMv - v'))x*Av'W, (58) 




with V k (r,) = lj 2 _ , k - = -(£ - l/6)C( V )R(r,). 

Note that these modes satisfy 



e -t|k|T e -j|k|r 

* k ^ T ) = ^ffT and Xk,.(r)=-i|k|— ==. (59) 



Now applying Picard's method to lowest order one obtains 

/3 u (r; oo) - J 9{r - r^e-^V^dn, (60) 

where we've introduced the Heaviside function 9. Then from Plancherel's theorem one obtains the following density of 
produced particles at time r 



JV(r;-oo)^ (e lf f f T a^R 2 ^. 



(61) 



To calculate the rate of particles per unit volume and unit time, one has to multiply this last quantity by a 3 and to take 
the derivative with respect the cosmic time and finally to divide by a 3 , the final result is 

l-|(a 3 iV(r; -oo)) S -L(£ - l/6) 2 R 2 (r), (62) 
a 6 at loir 



that coincides with formula (2.29) of 02111 . 

Using the formula R = 12H 2 + 6^r, one has approximately at late times R = — ^j- s'm(Mt). Then averaging over a 
period of oscillation one finally obtains 



^aWfr-oo))*^-!/*)*. (63) 

At this moment an important remark is in order. If one makes, at late times, the approximation a(t) = i 2 / 3 (that 
corresponds to the expansion law of a matter-dominated universe) one has t ~ rf, then R(rj) ~ r/~ 3 what means that, 
at late times, Vk(rj) ~ rj and thus one has an infinite production of particles. One can avoid this problem neglecting the 
effect of power-law expansion and taking as a scalar factor the function a(t) = [1 + — sin(Mt)], in that case one has 
t ~ r\ and thus, at late times, one has Vviji) ~ V 1 - 

To calculate the energy density per unit time one has to take into account that the energy density at a given time, in 
general, diverges for a no-conformally coupled fields (see for example B56lO . then one has to re-normalize this quantity. 
Moreover, the re-normalized energy density depends, in general, on the regularization method used, and in the final result 
it appears some vacuum polarization terms, that are very difficult to separate of those relative to the particle production 
(only in the conformally coupled case this is possible to do). 

For these reasons in order to "calculate" the energy density per unit we use the fact that the main contribution to the 
particle production come from modes with k ~ M/2. This can be deduced taking into account the oscillating behavior 
of the curvature at late times, effectively, inserting Vk = — 1/6) — sin(Mrj) in the beta Bogoliubov coefficient, 
one deduces that the oscillating character of the integral disappear when k — M /2. Then one might concludes that an 
approximation for the energy density per unit time, averaged over an oscillation, is given by 



|,(r;-oo) * f ^(^(r;-oo)) * 1/6) 2 . (64) 
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At this point is important to remember that Starobinsky proposed in l2(ill that the oscillating behavior of the scale factor 
can be though as a coherent oscillations of a massive field described by particles of mass M (scalarons), then with this 
point of view, gravitational particle production could be understood as a decay of scalarons (due to its rapid oscillations) 
into other particles like massless no-conformally coupled particles, massive conformally coupled ones, etc... 

In this way, one can calculates the rate at which the energy of scalarons is dissipated 



P(r) 



(65) 



where to obtain the energy density of the scalarons p(r) one must uses the "effective" Friedmann equation H 2 = 8 jj^ P , 

being m p the Planck mass (we say "effective" because in the Starobinsky model the energy density that appears in the 
equations is only due to the vacuum polarization, but from Starobinsky's viewpoint, at late time, one could thinks that 
scalarons drive the expansion of the universe). 



Now since H = ^ one has 



3M 3 



2 m p 



2 (£-l/6) 2 . (66) 



Finally at time t ~ T 1 the oscillations of the scalarons field are damped, the created particles thermalize, and 
the universe becomes radiation dominated. To calculate its temperature, one has to use the thermodynamical relation 
p = ^qN(T)T 4 , where N(T) is the effective number of relativistic degrees of freedom at temperature T, the effective 
Friedmann equation, and the fact that reheating ends when T ~ H (see for example [57]), then one obtains 

T th ~ jf^~ p ~ If - l/6\M 3 / 2 m-V 2 . (67) 



To end the section, we study massive conformally coupled particle production in the Starobinsky model. First at all 
note that, to apply formula dTTl ) one must assumes that lim t _j. +00 a(t) — lim t _ i ._ 00 a(t) ll49ll . then one doesn't have to 
work with the Starobinsky model because it doesn't satisfy this assumption. We will work with the an scale factor that 



at early times is a(t) = 1 and at late time is a(t) = [1 



3MI 



sin(Mt)]. We also assume that m <C M, then since 



= |kl : 



m 2 , at late times, one has Vk(?y) = w 2 k - 0^(77) = m 2 (l - a 2 (r])) = sin(M?7). 



Now if one chooses the family of solutions 



one obtains 



/3 k (r;-oo) 



Vk(r)') sin(w_ !k (?7 - v'^XKrivldij' , 



8(t - ^e-^-^W^dr]. 



(68) 



(69) 



Thus, since now we cannot apply Plancherel's theorem, to calculate the density of produced particles at time r, we 



make the change of variable |k| 
one gets 



N(t; -00) 



= x 2 M 2 and make the approximation xM \J x 2 M 2 — m 2 = x 2 M 2 — then 



M 3 



2^ 2 a 3 (r) 
M 3 



/M 



x2 \fik(T; — oo)\ 2 dx - 



M 2 



,/M 



|/3k(r; -oo)| 2 da; 



27r 2 a 3 (r) J 
Now we can apply Plancherel's theorem to obtain 

N(t; -00) = 



x 2 \Pk{T; — oo)\ 2 dx. 



1 



167ra 3 (r) J 



V£(r))dr), 



(70) 



(71) 
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and then, at late times, one has 



l|( a 3 Ar(r; - 00 ))^_^ sin2 (Mr), (72) 



where we have used the approximation a(r) = 1. Averaging over an oscillation one has 



^l^-oo))-^^. (73, 

To calculate the energy density per unit time, we use once again that the main contribution to the particle creation comes 
from modes with k ~ M/2, then 



d , , Mid,. , , . m 

5>< T! -°°> " Y^Tt^ S 36^W (?4) 



And finally, the rate at which the energy of scalarons is dissipated is 



^p(r;-oo)_ m 4 



r = at — - = (75) 



this is the result obtained by Starobinsky [20]. 



B. Particle production by strong electromagnetic fields: Schwinger's formula 

Here we deduce the well-known Schwinger formula 12411 for spin and spinless particles using the W.K.B. approximation, 
that is, we compute the probability that the vacuum state remains unchanged in the presence of a constant electric field 
using the semi-classical approach. 

First, we consider, in the Minkowski space-time, the Klein-Gordon field in a box of volume i 3 , coupled with an external 
uniform vector potential f (t) . The Klein-Gordon equation is equivalent to a Hamiltonian system, composed by an infinite 
number of harmonic oscillators with frequencies which depend on time. The mode equations are: 

Xk + Wk(*)Xk=0 with k€Z 3 , (76) 

where now = | 2 -|^ + ef (<)| 2 + m 2 (being e the electric charge). 

We assume that lim t _ i .-|- 00 Wk(i) — Wk,±, and we write the "in"-states as linear combinations of the "out"-states as 
follows 

Xm,k(t) = akX<mt,k(t) + AcXLt.kW, ( 77 ) 

then one has, A out ,k = a k A in , k + /^ALk- 

Let |nk) be the "in"-state that contains n particles in the k mode, and let |nk) be the "out"-state that contains n particles 
in the k mode. Then, it is easy to obtain the following relations fl]: 

|0 k ) = C k £ % |n k ), \^)=C % Y i [-—) ( 78 > 

n=0 Va k^ n=0 \ " k ' 

with|C-k| 2 = |C k | 2 = | ak |- 2 . 

From these relations, we deduce that the probability that a particle in the k mode is produced [2, 6;|58[], is 

K-KWI" = (1+ '^V +1 . 09) 

and the average number of produced particle in the k mode, is 

<0 k |a+ k <w )k |0k) = |/3 k | 2 . (80) 
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Example II.3. (Adiabatic approximation) From formulae ( 1291 ) and ( 157 1 ) one obtains the following formula for the number 
density of produced particles 

N diag {t ] -^) = -^—\E{t)\, (81) 

where a = e 2 is the fine structure constant and E(i) = — f (t) is the external electric field. However, their energy density 
diverges. To obtain a well-known defined quantity one has to re-normalize the electric charge, after this one has ( see 
details in il EHzIl/) 

p diag (t; -oo) = ,jf 2 4 |E(t)| 4 + — (|E(0| 2 ~ 2E(i) • E(t)) . (82) 

Moreover it is also possible to calculate, after charge re-normalization, the induced electric field, that is, the corrections 
to the external electric field due to vacuum fluctuations |6^/ 



and the effective Lagrangian density for the electric field (i6 

C eff (t) = -^|E(i)| 2 + J Q2 2 j Eft)| 4 + ^^|E(i)| 2 , (84) 

which generalizes the Euler-Heisenberg formula KI^l . for a time-dependent field. 

From these results obtained above one can deduce that the probability that the vacuum state remains unchanged, namely 

Pis 

p = n t+jbj? = exp (- £ iog ( x + |/?k|2 )) - exp f - e £ ^i^i 2 ") < 85 > 

kGZ 3 '^ k| \ kGZ 3 / \ kGZ 3 n=l / 

As an application of this result, we'll find Schwinger's formula for scalar particles. Consider the case f(t) = 
(0,0, /(*)), where 

-ET if t < -T 
f(t) = { Et if -T < t < T (86) 
ET if t> T, 

being E the electric field and with T ^> 1. We suppose for example eP > (The case eE < is analogous). 
The Schwinger formula gives the probability that the vacuum state remains unchanged. Then, using the notation 



Schwinger's formula for spinless particles is [2 



I °° f-l)" +1 \ 
P = exp ~^2J — 2 — exp(-nS) . (88) 



In order to deduce this formula, we compute the /? Bogoliubov coefficient using the relativistic tunneling effect | 
i.e., using the formulae given by the W.K.B. method in the complex plane Il62l464ll in a formal way. (we say in a formal 
way because our field is not analytic, and the formula ( |89l that we'll use is only mathematically justified for analytic 
fields). The result is JH 



\M 2 =\ ^[-ImJ^^l+m^ + ^+eEz) )dzj if \^\ < eET (g9) 

if | ^ | > eET, 
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where k± = ^(fcijfe), and 7 is a simple curve in the complex plane, containing the complex turning points 



L ± ^ k±J —— as interior points (Note that for l^rH > eET there isn't turning points l66ll . that's the reason why 

for that modes its beta-Bogoliubov coefficient vanishes). Now, it's easy to verify that 



exp (_ ^( fc2 i+ m2 > ) if \!^\ <eE T 



if I ^ I > eET, 



and therefore, the probability that the vacuum state remains unchanged is 

/ 00 (_-\\n+l \ / 00 (_-,\n+l \ 

P = exp - E E HH^" = exp -iV E M,— exp (-nS) , (91) 

in agreement with Schwinger's result. 

In the same way one obtains that the average number of produced pairs per unit volume and per unit time, is 

E 2 a ( nm 2 \ / \ 

8^ eXP ('^)- (92) 

To obtain Schwinger's formula for fermions, one has to use the Pauli Exclusion Principle to get the following relation 
between Bogoliubov's coefficients II651 167[ 16811 

|«k| 2 + |/3k| 2 = 1. (93) 

Now I /5k 1 2 is the probability that a particle is created in the k mode, because it's the same that the average number of 
produced particles in that mode, and consequently |«k| 2 is the probability that no particles are produced in that mode. 
Therefore, due to the existence of two states with spin 1/2, the probability that the vacuum state remains unchanged is 

P=U( 1 l^l 2 ) 2 = exp ( -2 E E ^xl 2 " ) = ex P ( ~ 2Af E \ exp(-nS) J . (94) 

Remark II.5. In that Section we've obtained Schwinger's formula in an easy but f ormal way. However if one wants a 
complete demonstration of the formula with all the details, one can look up at i l72l475l/ . Essentially, the idea is to make 

the change of variable y — yj ( 2 ^° 3 + eEt), then the Klein-Gordon equation behaves 

Xk + ^k(y)Xk = 0, (95) 



where, w k (y) = \J ' \y 2 - and = ^(kj_+ to 2 ). 



In this case, the function Xk-0-WKB(y, e = 1) = \/ 2zt(y) e ' ' JJk ^ dT , has the asymptotic behavior 



e iy*/4( y 2s-l/4-iA k /2 ^ 

Xk,o-,wKB{y; e = l) = <j e _ i!/V4(y2) _ 1/4+iAk/2 y ^ ^ (%> 



Note that o-W-FCB^i £ = 1) = ~*^k(y)Xk,0;W k b{u] e = 1) when y — > ±00, this means that the diagonalization 
metod and zero order adiabatic vacuum modes give the same value for /3k- 

On the other hand, a independent set of solutions of^5\ is given by the two following parabolic cylinder functions ^2 

u K1 (y) = exp f-~y 2 J M (- l -A k + i, i, '-y 2 ^ (97) 

u k ,2(y) = exp (~^ 2 ) V exp M (-~A k + 1 ~ ~y 2 ) , (98) 

where M is the Rummer's function. 
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We now define the mode solution 

Xk (y) = ^/2 jB -l e W8 e -^ k /4 2 -l/4-M k /2^ k(2/)) (99) 

with 

My) - T % + { ^ Ak) u k Av) + r (l r ^ k) ^ 2 Q/), (ioo) 



r(l-|^ k ) r(|-|A k ) 



ant/ 



Then, from the asymptotic behavior of the Rummer's function [47] we can see that 

Xk(y) XKo-,wKw(y; e = 1), y -t -oo, (102) 

and 

Xk(y) -> B k 1 e -/4 2 i-^ Xk . WH7 (j / . e = 1) + C k B k - 1 Xk,o ; wxw(y; e = 1), vv/zen y -> oo, (103) 

r _ rQ + iA k ) , r(| + i^ k ) 
k "ra-M k ) r(i-i^ k )- 

Then from this last formula we can deduce that, in both cases (diagonalization method and zero order adiabatic vacuum 
modes), the square of the (3-Bogoliubov coefficient is given by 



|/?k| 2 - \C k /B k \ 2 = e 2 ^ = exp + m 2 )j , (104) 



in agreement with ( I90P . 



C. Moving Mirrors 



Consider a massless scalar field <j> in the 2-dimensional Minkowski space-time interacting with a perfect reflecting 
moving mirror. Assume that the mirror trajectory follows an inertial prescribed trajectory x — git), that in light-like 
coordinates u = t — z and v = t + z, we write as v — V(u) or u = U (u). 

For a perfectly reflecting mirror the set of "in" and "out" mode functions is |p76l — TtsIi 



V 47r M 



V 47r M 



e -^(«)) 0(y - V(u)) 
e -iwu{v)^ o{u-U(y)) 



(105) 



— e 



U ^)6{v -V{u)) 
e~ iuV ^) 6(u-U(v)), 



(106) 



, R (resp. 0™ L ) represents particles with frequency oj coming from the right (resp. left) null past infinity domain 
(resp. i7£")> and </>° u jij. (resp. 0° u £) represents particles with frequency cj going to the right (resp. left) null future 

1 T+\ 



where 
Jr (" 

infinity domain J R (resp. )■ 
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It is a well-known fact that the average number of particles in the ui mode produced from the vacuum in the right hand 
side (rhs) of the mirror, is 178l 17911 



N(lo) = J du>\d 2 . 
Then our main objective is to calculate the beta Bogoliubov coefficient 

<^(<k*;$U)*, with u , u '>o, 

where the parenthesis in the right member denotes the usual product for scalar fields Hill , that is, 



(107) 



(108) 



& = * / *) dt 0, R (t, x)dx, 

J 9(t) 



(109) 



that doesn't depend on the space-like hyper-surface chosen. 

The best way to perform this calculation is to choose as hyper-surface the right null future infinity J^, then one has 



& = 2z / d«c3AC« 



due- l " u d u e- llJJ ' v{u) = —\Hr I due- tuju e~ tuj ' v{u) . (110) 

Z7T V ijJ ' 



Assuming that the mirror's velocity converges fast enough to a constant when |u| — > oo, and integrating by parts one 
gets 



8* 1 ',, = — - — Vlju' I du 



1 



V"( u ) -iuiu -iu'V(u) 



{co + uj'V'(u)y 



(111) 



For simplicity we assume that the mirror's acceleration is discontinuous at the point u = a, after another integration by 
parts one obtains 

& = -^ ^^ e "^ yW(y/,(0 " ) ~ V " {a+)) 



2vr 



du 



V"'(u) 3lu'(V"(u)) 2 



{lj + lj'V (it)) 3 (lj + lj'V'(u)) 4 



—iu}U„—iojV(u) 



(112) 



From this formula, assuming that the mirror's trajectory is asymptotically inertial, i.e., V'(u) > Vm 6 1 (see for 



example [79]), one concludes that 



R,R 



and lj 



B R ' R , 



are integrable functions in the domain [0, oo) 2 \ [0, 1] 



Finally, we're interested in the production of particles in the infrared domain, i.e., we want to calculate 
[0, l] 2 . We write the Bogoliubov coefficient as follows 



nR,R 



RR = 1_ jw_ 

Pu < u ' 2ir V u>' 



with B > 0. After an integration by parts one obtains 

1 vww 



aR,R 



2ttlj + Blj 1 



7 / du{V'{u) - B)e-^ u e^ v{ - u \ 



(113) 



(114) 



and thus, if the function | V'(u) — B\ is integrable in R for some B > 0, it can be deduced that 
are integrable functions in the domain [0, l] 2 . An example of this kind of trajectories is 



pR,R 



and lj 



nH,H 

" UJ.UJ 



Bu u < 

V(u) = { V(u) < u < u 

V(uq) + B(u — uq) u > uq. 



(115) 
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However if we one is only interested in the convergence of the function ui 
trajectories that satisfy 



nR,R 
uj.uj' 



in the domain [0, l] 2 , one only needs 



/0 roo 
du\V'(u) - Bi\ < oo and / du\V' (u) - B 2 \ < 
-co JO 



(116) 



for some no-negatives constants B\ and B%, (here it's important to remark that one of these constants can be zero, that is, 
it's not worth that the trajectory be asymptotically inertial). To prove this statement we write 



gR,R = ! 



due^ 1 ^^ - B ^ u ') u e ~ iu] '( u )- B i u ) _|_ / fi ue -i{w+BW)u e -iu' (V(u)-B 2 u) 

Jo 



(117) 



and we assume for simplicity that V(0) = 0. After an integration by parts one gets the formula 



1 ^ 
1 



B\ — B 2 



1 



2ttlo + B 2 ui' 



(u + B 1 lo')(uj + B 2 uj') 2ttuj + Biu' 



du(V'{u) - B 2 )e~ lu,u er lbJ ' vl - u \ 



du{V'{u) - B x )e 



— iuju^-iuj'V(u) 



(118) 



that proves the statement. 

In conclusion we've proved that for asymptotically inertial trajectories with continuous velocity the radiated energy is 
finite. However it is also possible an infinite production of particles with very low frequency (an infrared divergency). To 
remove this divergency one must assumes that the initial and final mirror's velocity is the same. 

Note that particle creation for partially transmitting mirrors is a bit different: At very high frequencies the mirror 
behaves transparent, and then there are not particle production independently of the mirror's trajectory. On the other hand, 
at very low frequencies the mirror behaves a perfect reflector and then the same kind of infrared problems as in the perfect 
reflector case remain. Consequently, if one is only interested in a finite radiated energy, one can considers trajectories with 
a continuous velocity V'(u) Vu € K, which fulfill the condition ( II 16b as, for instance, the no-asymptotically inertial 
trajectory: 



v = V(u) 



u if u < 

i(l - e~ ku ) if u > 0. 



(119) 



1. Simulating black body collapse 

Now we're interested in a trajectory that simulates a black body collapse [11, 80, 8l]], that is, with the following form 

{u if u < 

V(uq) + A(u — u ) if u>u Q , 

with A = e~ kua , where k is a frequency and kuo ^> 1. 

Note that this trajectory can be written under the following form, too 



v if v < 

u = U{v) = { -i ln(l - kv) if < v < v 
U(v ) + A' 1 ^ - w ) if v > v . 



(121) 
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Then, 



1 



1 



_ g — iwuog— i(/V(u ) 



1 /u/ 1 
2tt V w" 



l-A 



ds(l-s) lw/fc e" lw '/ fcs . (122) 
Assuming for simplicity to ~ A;, if one is interested in the domain of frequencies 1 <C w'/fc <C A , one arrives at 

(123) 



pR 7 R 

Puj.uj' 



2ni 



^ww' 27r V u> k J Q 



) / 'k g — iui' j 'ks 



To obtain an explicit expression to the second term on the rhs, we consider the domain 

k 

D = {zE C/Rez G [0, 1 - A],lmz G [-e,0], with — < e < 1} 
and, going through the same steps as in ll82lh one easily obtains 



2iriy/<jj(jj' 

and thus, using that |T (1 + iuj/k) | 2 - ™ /k 



ik 



j Ik 



— r(i + iw/fe), 

u' 1 



sinhf-Trcj/Zc) 
2 



(see B47I1 ) one gets, for a perfect reflecting mirror, that 

oR,R 2 ~ 1 f p 2™/A _ lA" 1 



(124) 



(125) 



(126) 



in the range 1 < oj'/k < A^ 1 . 



Remark II.6. From that formula, one deduces that the number of radiated particles in the U) mode diverges logarithmic- 
ally with uq — > oo. In this situation the physically relevant quantity is the number of created particles in the oj mode per 
unit time. This dimensionless quantity is finite and its value is given by rt&2l l&jl/ 



lim 

jio— s-oo Uq 



1 -N(lj) = — ( e 2 ™' k - l) \ 

In V ' 27T V J 



(127) 



Now we'll study what happens when the mirror is partially reflecting. First, we search for the co-moving coordinates 
(t, p), that is, the coordinates for which the mirror is at rest, r being i the proper time of the mirror, and we take p such 
that its trajectory is given by p = 0. Introducing the light-like coordinates (u, v) defined as 



u = r — p\ v = t + p, 



(128) 



we will calculate the mirror's trajectory in the coordinates (u,v). Along this trajectory, the length element obeys the 
identity B 



dr 2 = du 2 = dv 2 = V'{u)du 2 = U'{v)dv 2 . 
Then, an easy calculation yields the relations 

u if u < 



(129) 



v — u(u) — < 



|(1 - e r ku / 2 ) if < u < u 



(130) 



u{uq) + V A{u — Uq) if U > Uq, 



and, 



if v < 



u — v(v) = < |(1 — \/l — kv) if < v < v 
v(vq) + A~ x l 2 {v — vq) if v > vq. 



(131) 



20 



When the mirror is at rest, scattering is described by the S-matrix (see 

s(uj) r{uj)e~ 2iulL 
r(w)e 2lu)L s{uj) 



for more details) 



(132) 



where x = L is the position of the mirror. The S matrix is taken to be real in the temporal domain, causal, unitary, and 
the identity at high frequencies |86]. Correspondingly, the "in" modes in the coordinates (u, v) are |[88ll 



1 



' s(w) e - MC 6»(u -v)+ / / | (e~^ c + r{uj)e~ luU ) 6(v - u). 



(133) 



(134) 



On the other hand, the "in" modes in the coordinates (u, v), namely </> m , are defined in the right null past infinity JZ 



by 



and in the left null past infinity J L by 



-tujv . ±in 



0, 



(135) 



(136) 



Writing g"\(u, v) = g 1 ^ ^.{u(u) 1 v(v)) with k = R,L, and using that g™ w fe = .g" 1 ^, one obtains the following relation 



Cfc = dw'x(w')(ff", fc ; Cfc)ff™,fc. with k = R,L 



(137) 



with x( w ') the sing function. To calculate explicitly the "in" modes, I choose the coefficients r(w) = jj^- and s(uj) = 
with a > that correspond to the so-called Barton-Calogeracos model ll88l - l90ll . In this case, on the rhs of mirror 



one has [91] 



where 



i 



e -iM + ^ {u) . m L ( Uj v) = ^trans ^ 



(138) 



C2(«) = < 



1 -i7_ e -7ii(u) 



£7_ -iuV(u) . 

1 



fe^47r|aj| 



w+27 



e 



u < 
< u < uo 



(139) 



CJ+Z7 



_L_ r e -ta»V(ti) _ e -ta»V(ti D ) e -7(tt(«)— u(u ))j 



-^/47t|lj| V3"+i7 



27 e -»f e - 7 (g(u)-u(uo)) fl«(«o) dse i«( s+ i_i ii ( U0 )) 2 e . 



and 



u < 



1 



+ 



1 -»7 p-ew(u) 27 r#S(«) . / , 1 _ fc-/ ..\\2if n < < 



±L e -7»'(") J. , 1 e ,1J " 
fcy47r|w| 



K«)-«(«o)) + j 7v A4 e -7(«(«)-G(«o)) 



2 7=e - 7 («(«)-fi(«o)) fl«(«°) ds(s + 1 - |u(u )) 2 ^e^ 



(140) 



U > Uq 



21 



Note that in the case of perfect reflection, that is, when 7 — > 00 one has 



1 



3 — iwV(tt) . itran 



and when the mirror is transparent, i.e., when 7 — > one has 

0; 



Since we're interested in the particle production on the rhs of the mirror, we must now calculate 

%*;^,nT, and CMc;r;C,if u,,u'>0. 



oR,R — / lput *, ^in 



(141) 



(142) 



(143) 



In order to calculate this products we choose the right null infinity because here the "out" modes have a very easy 
form, then 



We start calculating /3%'*, = 2% J* R ducl)° u R d u fJ f R , with the result 



(144) 



1 



7 



27ryW/ w + ?7 



1_2 / ^-A-i/s^a-v^/fc 



_| 7 p —iw'/k / ^ xx iu)/k-l/2 e iu'x/k 

2irki^f 'louj' J a 



27 
& ./o 



e i:w'(s 2 +2sV^)/fe e -27s/fc 



(145) 



Now assuming once again uj ~ fc, provided that 1 <C 77 <C ^2- <C ^4 _1 , equation ( 11451 ) turns into equation ( 11231 ). 

2 

Consequently, we precisely obtain the same behavior as for a perfect reflecting mirror. However, in the case 1 -C ^ <C 
t- <C A -1 we observe that 



R,R 



2ixki\/u)Lo' 

12 _ 7T 



and using the identity |r (1/2 + iw/fc) | - cosh(Wfc) 

2 1 / /y \ - 



iui/k+l/2 

r(l/2 + iw/fe), 
(see B47H ). we conclude that 



1 



2ixku) vcj' 

Finally, a simple but rather cumbersome calculation yields in the first case 



ii.L 



and in the second one 



oil,!, 



2 1 



0. 



70 



(3V 



(146) 



(147) 



(148) 



(149) 



Then we can conclude that the number of produced particles in the mode uj is approximately given by |92, 93 



N(uj) = 



i 2 /k 



^_(Mk_ 1 ) 1 ,f _J_(Z\ 2 ( e ^/k + 
2irkuj' V I J ry 2 /k 2irkio\uj'J \ 

^-bxh/k) (e 2 ™l k - lY 1 + J- (e 2 ™' k + lY 1 Si^-Uh/k) (e 2 ™/ k - lT\ 
Trk \ ) 2ttuj V / 7rk V / 



(150) 
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because we're assumig k ~ u <C 7. 

That is, the number density of produced particles in the mode to by a partially transmitting moving mirror is finite when 
Mo —> 00, moreover when lu ~ k <C 7 the mirror radiates a thermal flux described by Bose-Einstein statistics. However 
when cj ~ fc ~ 7, maybe the contribution of the sector [7 2 /fc, 00) could be dominant and another kind of statistics 
(Fermi-Dirac) would be possible. This is a situation that deserves futher investigation. 



III. VACUUM FLUCTUATIONS 



The re-normalization of the two point function via adiabatic regularization is given with all the details, and the re- 
normalization of the stress tensor is also reviewed. 



A. Re-normalized two point function 

1. Massless conformally coupled field 

We start this lecture studying the simplest case: a massless conformally coupled scalar field in a Friedman-Robertson- 
Walker space-time. We'll calculate the re-normalized part of the two-point function (</> 2 (?7, x)) using the adiabatic regu- 
larization method. This is the simplest example, and it help us to understand all the details of the method. 

Consider the quantum scalar field 



<£(j7,x)= / d 3 k Ok^k(»7,x) -Fa^C^x) 

JR3 L 



(151) 



-i|k|f) 

then, using the same notation as in first lecture, one has Xk\H) = / > an d thus, the two-point function is given by 

(0 2 (^,x))^/ ^k^x^—L^ f°|k|d|k|. (152) 

In order to obtain the re-normalized value of the two-point function, we'll follow Bunch's method described in l34ll : 
First, one considers the adiabatic modes obtained in the WKB approximation, 

XadiM (r,) = ^=e- i f w * d *>, (153) 

up to order 2. To calculate these adiabatic modes one has to use equation (0, then in the conformally coupled case Wk is 
given by (see for details H0B) 



W _ uk _J4 + »^ (154) 



with cj 2 . = |k| + C(rf)m , and a simple calculation yields, 

lm 2 C" 5 m 4 (C") 2 

Once the adiabatic modes has been calculated, to obtain the re-normalized expression of (<f> (t), x)) one has to subtract 
from ( 11521 ). the adiabatic terms up to order two (only terms that contain, at most, two derivatives of the scalar factor) that 
appear in the expression Ir96ll97ll 

Ikl 2 

-d|k|, (156) 



47T 2 C(77) J 
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and finally to take the limit m — > 0, that is: 

#{ „, X)U _ „ ' [ r f M _ m , ik! + rr^ r tt, |k| _ rtt, ik| i , 157) 

m->o A-K z G(ri) [Jo V w k / 8 J w£ 32 J w k 

It's not difficult to show that the final result is 

which coincides, for the de Sitter phase, with formula (3.19) obtained in [98]. Note also that, in the case of an universe 
filled by radiation, and consequently R = 0, (4> 2 (ri, x)) ren = in the massless case independently of the coupling 
constant value. 



~c 



1 a" 



1 



48tt 2 a 3 288tt 2 



R. 



(158) 



Massless minimally coupled field 



In this Section we consider another simple example, a massless minimally coupled field in the flat chart of the de Sitter 
space-time, where the modes can also be calculated exactly, they are 



Xk(v) = (ak-0k(?7) + bkipl{r))) 



(159) 



where ipiiiv) I s gi ven by formula d43l and a k and 6k are some constants. 

In general, (cj) 2 (v, x)) = 2 Tr' i c(r)) Jo°° |k| 2 |Xk(f?)| 2 d|k|, has ultra-violet and infra-red divergencies. To avoid these last 
ones (see for details iMH), we consider a transition from the radiation dominated phase to the de Sitter one, described 
by the following scale factor: 



a(r]) 



2 - V/Va V < Vo 
Vo/V V > Vo, 



(160) 



with r/o = —1/H. 

The modes Xk f° r V < Vo, are given by e _i l k l''. Note that these modes correspond to the usual choice of the vacuum 
state for a massless field in the radiation phase, because in that phase the scalar curvature vanishes, and then Xk satisfy 
the equation x'k + |k| 2 Xk = 0, and consequently, the vacuum state is obtained in the same way as in the Minkowskian 
case, that is, from the modes e~ J l k l''. 

Matching at the point r\ — r/o the modes an their derivatives, one obtains 



1 



H 



H 2 



H 2 



re H 



2i|k| 



2|k| 2 ' 2|k| 2 " ~* ' 3H 

With these coefficients, for small values of |k|, in the de Sitter phase (rj > r/o), one has 

2 n H 2 V 2 



0(\k\ 2 /H 2 



IXk| 



1 

2|k| 



3Hrj 



(i 



C (|k| 2 /iJ 2 ) 



(161) 



(162) 



what shows that there is not infra-red divergencies. 

To analyze ultra-violet divergencies we calculate for large |k| 



Xk 



1 

2|k| 



1 



1 



H 2 



cos(2|k|(i?- 1 + v ) + 0(H 3 /\k\ 3 ) 



|k| 2 77 2 |k| 2 

what shows that the terms that give ultra-violet divergencies in (cf> 2 (i], x)) are 

V 2 H 2 



4tt 2 



t H f 1 



(163) 



(164) 
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Once we have separated the divergent terms, we calculate, up to order 2, the adiabatic terms that for 7] > rj o, are given 

by 

1 lm 2 C" 5 m\C) 2 

w *=^-^-s^r + T2^r- (165) 

From that, one sees that the divergent parts of ( 1156b are, in the de Sitter phase, given by 

Subtracting, for large frequencies, the divergent part of (1156b from (1164b , for example for |k| > \ri\~ 1 , one gets 

„2 tt2 roo / | k | 2 \ ff 2 r°° / i lid 2 \ 

lim ^TT / U k l - — d M = and lim / TTT - = °> ( 16? ) 

what shows that the ultra-violet divergencies are canceled. 

The problem now is that the subtracted adiabatic term J^ 1 J^-d|k| contains an infra-red divergency, be- 
cause it diverges when the mass approaches to zero. It is important to remark that this term can be written as 

"sf 7 So (£ _ \) ^r d \K that is ' il does not 

appear in the conformal coupled case. 
The solution of this infra-red divergency emerges from the following observation: The adiabatic approximation is based 
on modes of the form (11531 ). and it is clear that this form only has sense if the exact modes, namely Xk, are oscillating, 
that is, if Xk satisfy the equation ©, with 

n£(v) > o. (168) 

In our case this conditions means |k| > -\/2|?7 _1 | = \2He Ht , and thus, we must only subtract adiabatic modes well 
inside in the Hubble horizon at time t. Consequently, no infra-red divergencies appears. 

Here an important remark is in order: Our recipe to eliminate the infra-red divergency does not affect to the conservation 
of the renormalized stress tensor which for a FRW metric reduces to (pC 3 / 2 )' + p{C 3 / 2 )' — 0, where p is the energy 
density and p is the pressure, because the adiabatic regularization lies in subtracting adiabatic terms up to a given order 
and then the conservation equation is safisfied for each order, and more important, this subtraction can be performed mode 
by mode SHI. This means that if one denotes by p a d(4>k) and p a d{4>k) the adiabatic terms of the energy density and 
pressure for the adiabatic mode defined in equation ( 1153b . the conservation equation (p a d(</)fc)C 3 / 2 )' + p a d(4'k)(C 3 ^ 2 y = 
will be satisfied, and since the substraction is performed mode by mode, one can subtracts a given number of modes 
mantaining the conservation equation. 

Summarizing, the re-normalized quantity is given by 

1 1 r °° |2| M2 - ■ 1 ^° 



<0%x)) re „ = lim _ — - / |k| 2 | Xk (r?)| 2 d|k| - - / i-Hk| . (169) 
m-yo 2vr 2 C*(77) \J 2 J^2 /M W k J 

Now since, the ultraviolet divergencies are cancelled, and 

Um r ^Kd|k| = lim r ^S^d|k| = 0, (170) 

using ( 1163b we finally get 



2 tt2 rs/2/\n\ 

(0 2 (r?,x)) ren =\^\ |k| 2 |xk(^7)| 2 rf|k| 

2r/2 r°o / tt2 



/oo / tt2 \ 

r% | C ° S(2|k|(jff_1 + 7?) + ° (^ 3 /|k| 3 ) J |k|d|k|. (171) 

If we are interested in the late time behavior (\r)\H <C 1, i.e., Ht ^> 1), we can make the following approximation 

2o-2 r V2/M i r H h 2 r^/M i 

(^ (7?)X )) ren ^L^ | k | 2 | Xk(r? )| 2 d | k |c,_^ | k | d | k | + _y _ d | k |, (172) 
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where, in the first integral we have used the approximation |xk(?7)| 2 — 9r/ i ggjE ( e 1- (I162l i). and in the second one 
Finally, after integration, at late times we get 

(t 2 (r h x)) ren - + In 2 + ff*) S J^t, (173) 

that coincides with the early result obtained in 11261 1351 l36ll . 

Note also that formula (11681 ) justify the prescription given in J27L l99ll . where the authors assumes that only modes 
outside to the Hubble horizon at time t contribute to the value of the re-normalized two-point function. 

In the opposite case, that is, for a few Hubble times, (for example t = l/H), it is not difficult to show that 

(0 2 (ry,x)) ren ~0(if 2 ). 

Remark III.l. In inflationar y cos mology the re-normalization of the two point function is sometimes obtained in a 
different way (see for example klOl\l ). The modes Xkiv) are given by formula (|43p , and thus 

To avoid the infra-red divergency one assumes that the initial size of the universe is of the order H , then the Fourier 
expansion shows the modes has to be wave-length smaller than the Hubble horizon, that is, impose the cut-off |k| > H. 
And to avoid the ultra-violet divergency one has to subtract adiabatic modes that satisfy equation (|76(S[) , i.e., adiabatic 
modes well inside in the Hubble horizon, and consequently one only has to take into account the modes that leave the 
horizon, more precisely modes that satisfy H < |k| < y/2/\r}\, then from (|7<57[) and M70\ one has 

1 /V2/H / 1 \ 

^ )) = ^mL |k| ( 1+ w) d|k| - 

Finally note that the first term is the usual contribution from vacuum fluctuations in Minkowski space and must be 
eliminated by re-normalization, which give the following re-normalized two point function 

I rV2/\i)\ 1 H 2 1 H 3 t 

when Ht > 1. 



3. Massive case 



First, we study the minimally coupled case with m <C H in the de Sitter phase. (This situation appears when the 
inflation fields is in the slow -roll phase, and scalar field fluctuations described by the two-point function are very important 
in order to understand the self -reproducing universes in inflationary cosmology (see for example IT26ll27ll )). In last Section, 
we have seen that the re-normalized two-point function is given by the formula ( |169t without the limit rn — > 0, because 

the adiabatic modes that we have to subtract satisfy |k| > -^yl — — -j^- 

The calculation of the two-point in the massive case is more difficult than in the massless one, however, at late times, it 
is possible to approximate its behavior very well. To do this, note that the adiabatic regularization method guarantees that 

f°° |k| 2 ( Ixk^l 2 - d\k\, (174) 

is convergent [43]. To perform this integral one can choose mode solutions that correspond to the Bunch-Davies vacuum 

state, i.e., Xk(v) = yj^e~ l ^ + ^ Hl?\\k\-q) with v = — ffi — | — jffr- Then, at late times, using the asymptotic 
expansion for large arguments of the Hankel functions (formulae 9.2.8-9.2.10 of l47ll ) one can shows that the divergent 
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terms of (11741 exactly cancels, and since a easy calculation proves that the convergent ones are of the order 0(m 2 ), one 
can disregard its contribution to the two-point function. Thus, at late times, we have 



(0 2 (?7,x)) 7 . en 



2 tt2 r H 



|k| 2 |Xk(r?)| 2 d|k| 



2f7-2 r V2/W 



rfH_ 
2tt 2 



|k| 2 |XkW| 2 d|k| 



(175) 



H 



In order to avoid infra-red divergency, we can calculate the first integral assuming a phase transition to the radiation 
dominated universe to a de Sitter phase at time ?/o = —1/H. In fact, we consider the general mode solutions in the de 
Sitter phase 



Xkfa) = V^V 1/2 (a k H& (|k|r?) + & k ff«(|k|r 7 )) 
and match the modes and their temporal derivatives at rj = —1/H, to obtain 

If 



(176) 



2i 



1 /7r|k|?/o 



bk = „ 

k 2i V 2 

Then from the small-argument limit 



1 Tr\k\ Va 



2|k 



ff«(|k|,fe)-2rW'(|k|.») e^/* 



Hl 2 \\k\ Vo )^-HM(\k\ Vo )^-r(v) 



Mm 



and using that (see [47]) 



H 



(177) 



(178) 



(179) 



one arrives at the result |xk(f?)| 2 — gjc| (-^l 7 ?!) 1 2l '> an( l finally one obtains 



2H-2 pH 



rfH_ 

2tt 2 



|k| 2 |XkW| 2 d|k| 



Jl. 

18tt 2 



e 3H 



=fT when l/.ff < i < i7/m 2 
when t > H/m 2 . 



(180) 



with agrees with the first term of the right hand side of ( 1173b . 

The second one, can be done using the following approximation, valid of |k| > H, 



(181) 



Effectively, from (11771 ) one can easily obtains, in the range |k| > H, a k = 1 and 6 k — 0. Then, since at late time we have 
H^(\k\r)) = — ^r(i/) f-^j^ , inserting this expression in the second integral of ( 11751 ). one obtains, for a massive 



minimally coupled field with m <C H , 



2tt2 r V2/\v\ 



H 



|k| 2 |xk(r?)| 2 *- 



3H 4 



2jt^ - (F|7j|)sh5 



3H 4 
8m 2 7r 2 



1 — e 3" 



(182) 



because we have assumed m <C i?. 

Then, since the first integral in the right hand side of ( 11751 ) is smaller than the second one, depending on the value of 
m 2 t/H one has 



((^ 2 (77,x)) re 



3H 4 



1 — e 3H 



At when ll H < f < i?/m 2 



8?7J 2 7r 2 

which demonstrates, at late times, the formula (7) of ll26Tl . 



(183) 



3 if 



4-t when i ^> H/m 2 
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Remark III.2. In inflationary cosmology one can chooses Xk(^) given in formula ( 17871 ), then if inflation starts att^O 
and the initial size of the universe is of the order 1/H, after subtracting the adiabatic modes well inside in the Hubble 
horizon at time t, one only has to take into account the modes that leave the Hubble horizon. Finally one can uses the 
small- argument limit M79\ , which is equivalent to eliminate the Minkowskian vacuum fluctuations, to get 

2m 2 t " 

1 - e 3H ; 
because we are assuming m -C H. 



6tl 

Rm, 2 TT 2 



2is? - (H\r)\)f^ 



6n 

8m 2 7r 2 



Finally we calculate the re-normalized two-point function, for a massive conformally coupled field with m <C H in the 
de Sitter phase, given by 

„2 rj2 / roo i poo 1112 \ 

(0% x )) re „ = IL_^ |k| a |xk(»?)| a d|k|--jf ^d|k|J, (184) 



where Wk is given by formula (1 1 55b . (Note that, in the conformally coupled case we do not have to disregard any adiabatic 
mode.) 

At late times, with the same kind of argument used to disregard the terms that appear in equation dl74t , we can do the 
following approximation 

„2rr2 / MM" 1 , fA^ 1 

lllt^ —\\ 'I n II 1 1_ j 2 1 _ /„\|2jI1.I 1 



where A is some dimensionless constant of order 1. Actually, one can chooses A = 1, and the reasoning does not change, 
because \rj\~ 1 is large enough at late times. 

Since in this case there is not infra-red divergency, we can choose the modes corresponding to the Bunch-Davies 
vacuum state, that is, Xk(v) — \f^?e~ t ^ + ^H^ 1 (|k|^). Then for k < A/\r]\, a good approximation is given by 
Yk(f?) — — ?==, and since 

H 2 rf [ AH ~ X /„ . |k| 2 



4tt 2 



J (jk| - l -J-J d\k\ ~0(m 2 In (H/m)), (186) 



one can disregard this term (it is small compared with the other ones that are of order 0(H 2 )) and, at late times, we get 
the same result obtained in formula ( 11581 ) with R = 12H 2 , that is, 

(^ 2 (r y ,x)) re „ = (187) 
because the leading terms in ( 11851 ) are the same as in d!57| i. 



B. Re-normalized stress-tensor 



In this section we review some classic results about the re-normalization of the stress-tensor in FRW cosmologies. 
The vacuum energy density p is given by 



p = (T tt ) = (4tt 2 C 2 )- ] 



d|k||k| 2 (|x k | 2 +^|xk| 2 ) + 3(^-1/6) 



£(XkXk+XkXk)-^ 2 |Xk| 2 



,(188) 



where D = C'/C. 

To obtain the re-normalized value, one can uses the adiabatic regularization which consist in subtracting adiabatic 
modes up to order four. Following ll56ll one has to subtract the following divergent terms: 



(4^,-., dim ^, aizwgl r m K _ »« - y r^K 

167T 2 C^ J UJ k 167T 2 C Jr. UJ? 



(189) 
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and 



and the finite terms: 



(£ l l^ 2 .(j2D"D-MD a -27D 4 ) [°° d\k\^-, (190) 



128tt 2 C* 2 v 7n w2 



^ 1 ^Vd-V-^ 4 



384tt 2 C 28807r 2 C 2 \2 4 8 

+ - ^' 2 - 3D 4 ) + ^^(18^ 2 + 9D 4 ). (191) 

The simplest example is when one considers a massless minimally coupled scalar field, in that case the modes are 
Xk(r)) = then one easily obtains IT34T, [Tol Hoi 



/2|k| 



'-■« f ->~-SS>5'(i D " J '-I D "-!^l- 



In the flat chart of the de-Sitter space-time choosing the modes 



Xk(r?)=cJ^F( 2 >(|k|70, (193) 



with v = y | — — 12£ and C = e *('?' + 1), Bunch and Davies obtained in lll04ll using point-splitting regularization 
(which for scalar massive field is equivalent to adiabatic regularization |[56lll05lll06ll ) 



Pvac = (f tt )ren = — ^ ^ K + ( 12 ^ ~ 2 )# 2 ] 

-m 2 (12£ - 2)H 2 - ^m 2 H 2 - i(12£ - 2) 2 ff 4 + ^- J , (194) 



where tp denotes the digamma function. (In exact agreement with the previous resolt obtained in 1110711 using Schwinger- 
DeWitt regularization procedure ll76lo 

This result holds for all values of m and £ except in the massless minimally coupled case, because when m and £ are 
close to zero one has lll08ll 

eh™ 2 ^ 2 + ^ ~ 2)H ^ ^ (I - ") * i2S?i^r' (195) 

and thus the limit in ( 11941 ) gives different answers depending the way that m and £ approach to the origin. In order to 
calculate p vac in the massless minimally cou pled case using adiabatic regularization on has to consider the £ = 0, use the 
modes given in d43l and finally take m — > ll34l,l56tl . gettin 

29H 4 
P - = "960^' 

which can be also obtained taking £ = in ( 1194t and afterwards m — > 0. 

Remark III.3. It's well-known that the trace of the stress tensor corresponding to a massless conformally coupled field is 
zero, however after regularization its value is not necessarily zero. It was showed in [5^1 that the trace anomaly provided 
by adiabatic regularization, in the massless conformally coupled case is given by 

T vac EE (TS)ren = qqq^ i D '" ~ D ' D ^- ( 197 ) 

Then from the relation T vac = p vac — 3p vac one gets the value of the vacuum pressure p vac = (T xx ) ren . In fact, if one 
knows one of these values the others come from the trace anomaly relation T vac — p vac — 3p vac and the conservation 
equation (p vac C 3/2 ) + p vac (C 3/2 ) =0. 



29 



IV. AVOIDANCE OF COSMOLOGICAL SINGULARITIES 



The vacuum quantum effects due to a massless conformally couped field are taken into account in order to avoid 
classical cosmological singularities. 



A. Review of classical and quantum cosmology 

Classical cosmology: We will use the following notation: k 2 = lQnG = 16ir/m 2 , being G Newton's constant, p 
energy density, p pressure, lu a dimensionless parameter and H = a /a the Hubble parameter, being a the scale factor and 
the dot denotes denotes the derivative with respect to the cosmic time t. Then the Friedmann equation and conservation 
equation, for a fiat Friedmann-Robertson-Walker cosmology, can be written respectively as: 

H 2 = ^r P , p=-3H(p + p), (198) 
6 

and the equation of state for a barotropic perfect fluid, that we will consider in the paper, has the form p = top. 

With the derivative of the Friedmann equation, and the other two equation one easily obtains the "acceleration" equation 

H = ~(l + «)p. (199) 
Combining ( 11981 > and ( 11991 >, one can deletes p, and obtains the equation H = — 1(1 + uj)H 2 , then integrating one gets 

H{t) = 3(rby f^? (200) 

where t s = to — 3Hq ^ 1+lj ^ , being Hq = H(to) the initial condition. 

From the definition of the Hubble parameter, the following behavior for the scale factor is obtained 



( ^— ^ ) ■ (201) 



and from the Friedmann equation, one has 



Pit) = 3«»(l + o,)» ft-*.) 2 ' (202) 

The following remark is in order: if one assumes Hq > 0, then for u> > — 1 one has t s < to, that is, the singularity is at 
early times (Big Bang singularity), on the other hands, for to < — 1 one has t s > t , that is, the singularity is at late times 
(Big Rip singularity) ill ldll . 



Quantum Effects: Using the same notation as i39ll . for a massless field conformally coupled with gravity, one has the 
following expression for the trace anomaly 

T vac = 6a(H + 12H 2 H + 7HH + AH 2 ) - 12fi(H 4 + H 2 H) (203) 

with (see for example [ 109]) 

01 = 288^ (iV ° + 6Nl/2 + 12Nl) ' 13 = 288^ {N ° + UNl/2 + 62Nl) (2 ° 4) 

where A^o is the number of scalar fields, N 1 / 2 is the number of four components neutrinos and N\ is the number of 
electromagnetic fields. 

Remark IV. 1. The constants a and j3 come from the regularization process. For example dimensional regularization gives 
formula tO04\, and point-splitting gives (see [39]) a = 28 ^ 2 (N +3N 1/2 -18Ni), (3 = 2S8 ^ 2 (jV +ll/2iV 1/2 +62iV 1 ). 
Then, since the method of regularization influences these values, and it is uncertain what fields are present in our universe, 
one can considers all values of both parameters. 
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Remark IV.2. The relation between the notation of this section and the one of [21 ] that we used in section (II.A.3) is 

As explained in Remark III. 1 , to obtain the vacuum energy density one can uses the trace anomaly T vac = p vac — 3pvac, 
and the conservation equation. The result is 

Pvac = 6a(3H 2 H + HH- ^H 2 ) - 3(3H 4 , (205) 

which coincides with eq. ( 11921 > if one only considers scalar fields. Then taking into account this vacuum energy density, 
the modified Friedmann equation behaves 

K 2 

H 2 = —( P + p V ac). (206) 


With the derivative of this last equation, the conservation equation and the trace anomaly, one obtains the modified 
acceleration equation 

H = -^j" ((I + Uj)p + Pvac + ^{Pvac - T vac ) \ . (207) 

From both equations, one can deletes p, and one obtains the following third order differential equation 

4 • 1, „ , ,6 

n 2 Pvac 3 

that in terms of the Hubble parameter is given by 



- —H - p vac - -(p vac - T vac ) = (1 + lu)—H 2 - (1 + Lj)p vac , (208) 



-iff -(1+cj)4# 2 - 3(3{u + l)i? 4 + (18a(w + 1) - AP)H 2 H 

K K 



6a(u + 2)HH + 3a{uj + 3)H 2 + 2aH = 0. (209) 



B. a = 



In this section we consider the simplest case (a = 0), and we will see that, quantum effects don't avoid the singularities. 
Equation (fT2l reduces to the following first order differential equation 

H = ±- K \ , JL ■ (210) 



4 V ^ + PH 2 

1 . First, we consider the case a = and /3 > 0. Integrating equation ( fT3l one obtains 



- 1/H(t) + ] J^f arctan U ^-H(t)j = - |(1 + u){t t s ) + sj^- arctan (J ^H Q j (211) 

For u> < — 1, when t —> — oo one has H —> (that is quantum effects are small at early times), on the other 
hands, when t — > t s = t s — 3 ^ 2 +UJ ^ \j f tt — arctan (j ^piJ ^ ^ one has H — > oo, that is, the quantum 
effects don't avoid the Big Rip singularity that appears at t s . 

For u; > — 1, when t — > oo one has i? — > (that is quantum effects are small at late times), on the other 
hands, when t — > t s = t s — 3 ^ 2 +UJ ^ \j — arctan f^J ^^Ho^j one has H — > oo, that is, the quantum 
effects don't avoid the Big Bang singularity that appears at i s . 
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2. Now we consider the case a = and B < 0. The solution of equation ( fT3l is given by 



l/H(t) 



2H 4 



■In 



H(t)-H+ H + H^ 



H a -H+ H{t) + H^ 



= --(l+w)(t-t,), 



(212) 



with i/a 



• In the case w < — 1, for H e (0,iJ + /\/2) when i — > — oo on has H — > 0, and when — > H + /y/2 



at i. 



3(lH 



2H. 



In 



H+/V2-H+ Hq+Hj, 



one has 7J(i s 



-oo, that is, these 



H -H + H + /y/2+H + 

solutions are singular (the scalar curvature, R = Q(2H 2 + H), diverges). For Hq e (H+/y/2,H+) when 
t —> — oo on has — > however at t s H diverges. And finally for Hq G (H + ,oo) there is a singularity at 
finite time. Effectively , when t — > — oo on has H — » and when < — > t s on has iJ — > oo. 

In the opposite case ui > — 1, for i?o € (0, H + / v2) when Mooon has _ff — >• 0, and when — > Jf + / v2 



at 



3(l+w) 



H+/V2-H+ g +g + 



ffo-ff+ - H+ /V2+H + J ° ne h£lS ^) = +°°' * at iS ' theSe 

solutions are singular. For Hq € (H+ / y/2, H + ) when t —> oo on has H — s- however at i s H diverges. 
And finally for H £ (H + , oo) there is a singularity at finite time. Effectively, when t — > oo on has H — > H + , 
and when t — > t. on has — > oo. 



C. Empty universe 

Another simple case is the restriction to the invariant manifold p(t) = 0. In that case, one only needs the modified 
Friedmann equation, that is, the following second order differential equation 

1 K 2 8 

H 2 = n 2 a{"iH 2 H + HH - -H 2 ) ^-H 4 (213) 

Remark IV.3. The solutions with H > and those with H < decouple. To see this, we perform the change of variable 
Z = H/H to make the system no-singular at H = 0, then at H = the system behaves Z — Z 2 j% this means that the 
solutions can '? cross the axis H = 0. 

Equation ( |213t is an autonomous second order differential equation, then since solutions are invariant under time 
translations, the general solution is a one-parameter family of solutions. Taking this into account, we will prove that 
there is a one-parameter family of singular solutions. First we look for a particular singular solution, with the following 
behavior Hit) = j^t- near the singularity. Inserting this expression in (12131 ), and retaining only the leading singular 



terms one obtains C± — ^ ( — 1 ± y 1 + -§^J ■ Here it is clear that we have to impose the condition ^ > — 1. In terms 

of the scale factor one has ait) = clq ) ■ Then, for ^ > 0, the solution with C+ has a singularity of the type 

a(t s ) = 0, and the other one of the type a(t s ) = oo. However, when — 1 < <p- < 0, for both values of C, the solution 
satisfy a(t s ) — 0. 

Now we can prove that there is a one-parameter family of singular solutions whose leading term is H(t) = -r=h- To 

do this, we first transform the differential equation (I213l l in a first order one performing the change u(H) = Hit), then 
the equation becomes 

1 • k 2 8 

H 2 = K 2 a(3H 2 u + Huu' - -H 2 ) ^"# 4 > ( 214 ) 

where u'(H) = du/dH. Using the new variables H(t) = has the form u = —-^-H 2 , then if the equation is 
linearized about this point one obtains, for H — > ±oo, the following general solution 

u Hnearized = -^rH 2 + KlHl 1 - 30 ^ - when C ± ± 1/3 (215) 

G± K z a\&L± — 1) 
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and 

uunea^zed = -^~H 2 + Jjf - _L h | JT|, when C± = 1/3, (216) 
U± 3k z 

where K is an arbitrary constant. 

It's clear that we have to impose C± > —1/3. Then when — 1 < j- < or > 15, one has C± > —1/3, and then 
for both values one has a one-parameter family of solutions. When < ^ < 15, for C + , one has a one-parameter family, 
however for C_ one has to choose K = 0, that is, one only has a particular solution. 

Now we can perform the qualitative analysis. Firstly, note that when f3 < there exist two de Sitter solutions H± = 
±^J Making the change of variable p = \/\H\ (see i4lH42"lo . one obtains 

j t (f/2 + V(p)) = ~3ep 2 p 2 (217) 
where e = sing(H) and V(p) — — 4 ^ 2a (l + ^~P^j ■ in m the phase-space one has 

{l = -3ep 2 y y -V'(p). (218) 

We only consider the domain H > 0, because solutions with H > and H < decouple. The point p + = y/H+ is 
an extremum of the potential V, then, linearizing the system (1211 ). one obtains that for a < the critical point (p + , 0) is 
a saddle point, and for a > is a node stable. From the form of the potential, and taking into account that the system is 
dissipative in H > 0, the following results in the phase-space (p, y), with H > 0, are obtained: 

1. Case a > 0, (i > 

We have V < in (0, oo) and V(0) = 0. The (0, 0) is an unstable critical point. The solutions are singular at early 
and late times (p — > oo). Only a solution is not singular at late times, it is the trajectory that arrives at p — with 
zero energy (it arrives at the point (0, 0)), and only one if not singular at early times, it starts from p = with zero 
energy (it starts at (0, 0)). 

2. Case a < 0, (i > 

Now, V > in (0, oo) and V(0) = 0. The (0, 0) is an stable critical point, and solutions are only singular at early 
times. At late times they approach to the stable critical point. 

3. Case a > 0, < 

In this case, the system has two critical points. (0, 0) is an unstable critical point, and (p + , 0) is stable. Solutions 
are only singular at early times. At late times they oscillate and shrink around to the stable point, that is, (p + , 0) is a 
global attractor. Moreover, there is a solution that ends at (0, 0), and only a no singular solution that starts at (0, 0) 
(starts with zero energy) and ends at (p + , 0). 

4. Case a < 0, j3 < 

This is the Starobinski model lEoll . The system has two critical points. (0, 0) is an stable critical point, and (p + ,0) 
is a saddle point. There are solutions that don't cross the axe p = p + , these solutions are singular at early and late 
times, they correspond to the trajectories that can't pass the top of the potential. There are other solutions that cross 
twice the axe p — p + , they are also singular at early and late times, these trajectories pass the top of the potential 
bounce at p = and pass once again the top of the potential. There are solutions that cross once the axe p = p + , 
these solutions are singular at early times, however at late times the solutions spiral and shrink to the origin, these 
solutions pass the top of the potential once and then bounce some times in p = 0, shrinking to p = 0. These last 
solutions has the asymptotic behavior described in the beginning of section 2.1.3. Finally, there is only two unstable 
no-singular solution, one goes from (p + , 0) to (0, 0), and the other one is the de Sitter solution (p + , 0). 

Remark IV.4. Note that equation \16\ remains the same with the change H(t) —H(—t), this means that solutions 
with H < are the time reversal of the studied above. 
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Remark IV.5. For ^ < — 1 the values of C± are complexes, this is due to the fact that the system cannot go to (or come 
from) p — > oo monotonically, because the dissipation effect is not large enough compared with the potential force n4l\l . 

Remark IV.6. The case ui = — 1 (p = constant are invariant manifolds), is equivalent to the case of an empty universe 
with a cosmological constant. This case was studied with great detail in mt!/ . 



D. The general case 



The best way to study the general case is to consider the system 



H = Y 

Y = 1 ± w (2H 2 / K 2 -p/3-6aH 2 Y + aY 2 + f3H 4 ) (219) 
p = -3Hp{l + u). 



When (3 < 0, the system has two critical points (H±, 0, 0) whit H± = ±y — The semi-plane p = 0, H > (resp. 

p = 0, H < 0) is an attractor (resp. a "repeller") when ui > —1, and the roles are interchanged when u> < —1. 

What is important is to stress that in the case a < there isn't bouncing solutions, because at bouncing time, namely 
tb, one has —p(tb)/3 + aY 2 (tb) = 0, what means p(tb) = 0, but as we have seen in last section, p — is an invariant 
manifold where trajectories with H > and those with H < decouple. For this reason in the case a < there isn't 
stable no-singular trajectories, the only unstable no-singular solutions are the ones that appear in the Starobinski model. 

From this last paragraph one can concludes that, to obtain no-singular solutions, the interesting case is a > 0. In fact, 
the interesting one is a > and (3 < 0. 

In that case we can use the dimensionless variables t — H + t, H = H/H + Y = Y/H+ and p = ^Qr, then the system 
becomes 

H 1 = Y_ 

Y' = 2±H (-PH 2 f3p -6aH 2 Y + aY 2 + PH A ) (220) 
p' = -3Hp{l + uj), 

where ' denotes the derivative with respect the time t. 

In these variables the critical points are (±1,0,0). The linearized system at (1,0,0) has eigenvalues X± = 

—3/2 ^1 ± yj 1 + and A3 = —3(1 + uS). Since ReA± < and p = is an invariant manifold, all solutions in 

that semi-plane with H > go asymptotically towards this critical point. The eigenvector V3 = (1, —3(1 + u), 18u(l + 
u>)a/f3 — 2) corresponds to the eigenvalue A3, then for u < — 1, there is a solution that escapes to the de Sitter expanding 
phase following the direction of the vector V3. On the other hands, when u > — 1 the critical point is an attractor. At the 

other critical point the eigenvalues are A± = —3/2 ^1 ± \J— 1 + and A3 = 3(1 + oj). Since ReA± > and p = 

is an invariant manifold, all solutions in that semi-plane with H < escape from this critical point. The eigenvector 
i?3 = (1, 3(1 — 18cj(1 +uj)a/ (3+2) corresponds to the eigenvalue A3, then for u) < —1, there is a particular solution 
that goes asymptotically towards the de Sitter contracting phase following the direction of the vector V3. On the other 
hands, when cj > — 1 the critical point is a repeller. 

Now we look for singular solutions of the form H = C/(t — i s ), t — > tf. Inserting this value of the Hubble parameter 
in the conservation equation one obtains p(t) = p Q \i — t s \~ 3C ( 1+u \ where p has to be a positive parameter. Now , when 
uj < — 1, inserting the Hubble parameter and the energy density in the modified Friedmann equation one gets, once again, 

the values of C obtained in section 4.3, that is, C± = ^ ^— 1 ± + , because in that case the energy goes to zero 

when t — > if , and then, one obtains the same kind of result that in the case p = 0, the only difference is that now one 
has a two-parameter family of singular solutions (po is a free parameter, and the general solution of the system d23l is a 
two-parameter family due to the time invariance under translations). 
The case lj > — 1 is more involved. We summarize the results: 

1. When -1< £- < 
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• If C- < grrq^) . there are two two-parameter families with H = C± /(t — t s ) and po free parameter. 

• If C_ = 3 (i\ u s , there is a two-parameter family with H = C + /(t — i s ) and p free parameter, and a 
one-parameter family with H = C-/(t — i s ) and p~o = 0. 

• If C+ < 3 ( 1 + tJ ) < C_ , there is a two-parameter family with H = C+ / (i — t s ) and p~o free parameter, a 
one-parameter family with H = CL/(t — i s ) and p~o = 0, and a one-parameter family with H = 3 ( 1+ ^)( t _ t s 

ana p - - g ( 1+w )2 ^9[i+SP + T 3(1+^1 ~ Tj- 

• If 3 ( 1 ^_ tJ ) < C+, there are two one-parameter families with if = C± /it — t s ) and p~o = 0. 

2. When -1 = #■ 

• If 1 < 37j^j> there is a two-parameter family with H = l/(t — £ s ) and p~o free parameter. 

• If 1 > gnq^) » there is a one-parameter family with H = 1/ (i — t s ) and p = 0. 

3. When — 1 > there aren't singular solutions of the form H = C / (i — i s ). 

Remark IV.7. This result can be obtained in an equivalent way inserting in equation \209\ the function H = C/(t — t s ). 
Then, retaining the leading singular term, one obtains the values C± and gprp^j - Finally, transforming the differential 
equation in a second order one in the same way as we have done in Section IV, and linearizing around the singular 
behaviors obtained above, one can see the form and the number of parameters that depend the different families of 
singular solutions. 

To understand this summary we perform the change of variable p = y/\H\. Then the modified Friedmann equation 
becomes l42ll 

j t (p 2 /2 + V(p)) = -3epV + + U )p, (221) 

where V(p) = (l + ^p 4 ) - and e = sgn{H). 

The case cj < —1 is clear. Since p — >• at t — t s , one essentially obtains the same results as Section IV. However, for 
lu > — 1, on the right hand side of equation ( 12211 i one term is dissipative and the other one is anti-dissipative, moreover, 
in this case both terms diverge at i — i s . Then if one looks for singular solutions of the form H = C/(t — t s ), the first 
term in the right hand side of (l24l has to be dominant. And since this term is of the order l/(t — i s ) 4 , and the other one 
is of the order l/\t — i s | 3C, ( 1+ "), they will appear all the situations described above. 



It is also interesting to understand the form of the potential V (its picture appears in figure 3 of ref. 04211 ). It only has 
a zero at the point po = (3/2) 1 / 4 ^1 + yj 1 + | pj , and two critical points at p± = ( 1± v'J~ 4 £^ (p_ < p + y Then 
for p > 1/4 there aren't critical points, and the potential is strictly increasing from — oo to oo. For p < 1/4, the potential 
satisfy V(0) = — oo , V(oo) = oo and has a relative maximum at p_ and a relative minimum at p + (a hollow). For very 
small values of p at p_ one has H 2 = p, that is, the system is nearly to the Friedmann phase, and at p + one has \H\ = 1, 
that is, the system is near to the de Sitter phase. 

Next step is to find solutions that approximate to the Friedmann one when |t| — > oo (see |[54ll for the radiation case, 
i.e., cj = 1 /3). To do that, we consider equation (1209b in the dimensionless variables introduced above, and we reduce the 
order performing the change of variable u(y) — H(t) where y = H. The obtained equation is: 

2/3u + 3(1 + oj)/3{y 2 - y 4 ) + (18a(l + oj) - ^)y 2 u + 6a(2 + u)yuu + 3a(3 + u>)u 2 + 2a(ilu 2 + ii 2 u) = 0,(222) 

where now u = du/dy. Since the Friedmann solution in these variables is uf = — 1(1 + u)y 2 , the linearized equation 
about this point (u lmearlzed = u F + h) is 



2y-'-—h + - " y^ + Ay- 2 )h + B = 0, (223) 
1 + uj \9a{l+ui) z 



where A and B are some constants depending on the parameters a, j3 and u. 
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The idea to solve this equation is to take into account that for large values of |t| (small values of y), one can disregard 
the term Ay~ 2 . The homogeneous equation is solved performing the change h = \y\~ z, then one obtains: 

2 + q nl ^ z = °' (224) 



that we solve using the WKB approximation (see 115 511 page 276). Consequently the homogeneous equation has the two 
independent solutions 



hhomogeneous,±(y) = y 1/{ ~ 1+uj) exp (± ^_ y^-^/a^) 



(225) 



A particular solution is obtained using power series. It leading term is 



9a(l + u)) 2 

hparticular(y) = Jg By*. (226) 

Then the general solution of the linearized equation is approximately 

^Linearized Up "F ^■^'homogeneous, ± (y) ~t~ ^particularly) f^)r ^ y > 0, (227) 

where If is an arbitrary parameter, that is, we have proved that there is a one-parameter family of solutions that approx- 
imate to the Friedmann one for large values of \t\. 

Once we have seen these preliminary results, we can describe qualitatively the behavior of the no-singular solutions. 
We start with the case w < — 1. We have seen that there is a one-parameter family of solutions that at early times are in the 
expanding Friedmann phase, and we have to look for no-singular solutions that match, at late time, with that family. The 
only no-singular solutions at early times are: a one-parameter family that approaches asymptotically to the contracting 
Friedmann phase, and a particular solution that goes asymptotically towards the contracting de Sitter phase, following 
the direction (1, 3(1 + to), — 18u;(l + uj)a/f3 + 2). From the system d24l i it is easier to understand the dynamics. First, 
at early times the system is at the point p_ with H > 0. Then it leaves this expanding Friedmann state and rolls down 
either to the right or to the left. In the former case, the universe approaches to an expanding de Sitter phase (the relative 
minimum p + ). However since p is an increasing function with time, the critical points will disappear and the potential 
will be an increasing function with p, this means that the universe rolls down to p = 0, that is, it bounces and enters in a 
decreasing phase H < 0. Then it can arrive asymptotically at the points p_ oip + , (the no-singular solutions at late time), 
or it bounces many times in order to have enough energy in H < 0, to arrive at p — oo (singular solution). This last 
behavior can be easily understood, if one takes first into account that for H > (resp. H < 0) the system is dissipative 
(resp. anti-dissipative), and second that the energy of the system changes its sing when it bounces (see equation (12) of 

El). 

Finally, from the behavior of the no-singular solutions at late time, one can deduce that one has to very fine tune the 
initial conditions and the parameters a and /3 in order to obtain no-singular solutions that match these late time no-singular 
behaviors with the expanding Friedmann stage at early times, because these families of solutions aren't general solutions 
(a two-parameter family). 

On the other hand when cj > —1, we also have a one-parameter family of solutions that at late times are in the 
expanding Friedmann phase (in terms of the variable p this corresponds to the point p_ and H > 0), and we have to 
look for no-singular solutions that match, at early time, with that family. The only no-singular solutions at early times 
are: a one-parameter family that leaves the contracting Friedmann phase, in terms of the variable p, this means, that the 
system leave the relative maximum p_ with H < 0, and rolls down to the right or to the left, but in all cases, since the 
energy density is an increasing function of time in this region, the system goes to p = 0, i.e., it bounce and starts an 
expanding phase. The other no-singular solution at early times is a two-parameter solution (a general solution) that leaves 
the contracting de Sitter phase, in terms of the variable p, this means, that the system stars at p+ with H < 0, and then 
due to the anti-dissipation (at early times in H < the energy density is very small and the dominant term is the first one 
on the right hand side of equation d22U ) the system is released from the hollow and rolls down towards the region H > 0. 

From these no-singular early time behavior one can conclude that, in order to match these early time no-singular 
solutions with the expanding Friedmann phase at late time, one has to fine tune the initial condition. And depending 
on the values of the parameters a and (3 we will obtain different kinds of connections. For example, in lf55ll . different 
numerical calculations have been done in the radiation case, and they show the different connections in terms of both 
parameters. 
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